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Abstract 

In 1935 J.G. van der Corput introduced a sequence which has excellent uniform 
distribution properties modulo 1. This sequence is based on a very simple digital 
construction scheme with respect to the binary digit expansion. Nowadays the van 
der Corput sequence, as it was named later, is the prototype of many uniformly 
distributed sequences, also in the multi-dimensional case. Such sequences are 
required as sample nodes in quasi-Monte Carlo algorithms, which are deterministic 
variants of Monte Carlo rules for numerical integration. Since its introduction 
many people have studied the van der Corput sequence and generalizations thereof. 

This led to a huge number of results. 

On the occasion of the 125th birthday of J.G. van der Corput we survey many 
interesting results on van der Corput sequences and their generalizations. In this 
way we move from van der Corput’s ideas to the most modern constructions of 
sequences for quasi-Monte Carlo rules, such as, e.g., generalized Halton sequences 
or Niederreiter’s (t, s)-sequences. 
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1 Introduction 

If one is given the task to sequentially distribute points in the unit interval [0,1) such 
that for every non-negative integer N the initial N elements of the resulting sequence 
fill the interval very “uniformly”, then one might after some thinking come to an answer 
as illustrated in Figure [l] 

0 : 

1 : 

2 : 

3 : 

4 : 

5 : 

6 : 

7 : 

8 : 

9 : 

10 : 

11 : 



Figure 1: Intuitive construction of a sequence whose initial elements fill the unit interval 
very uniformly. 


That is, one might discover the sequence 

1131537195 13 
°’2’4’4’8’8’8’8’16’16’16’16’"" W 

A similar task is to distribute points from a sequence very uniformly on a circle. In 
this case, a possible idea may be to choose two antipodal points on the circle. Then 
rotate the whole picture by 7t/2 to obtain two further points. Next rotate the picture 
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by 7r/4 to obtain 4 new points. Next rotate by 7t/8, and so on (see Figure |2j) . In this 
way one constructs a sequence on the torus which directly corresponds to the findings 
in Figure [T| when the unit interval is coiled up into a circle. 



Figure 2: Intuitive construction of a sequence whose initial elements are uniformly 
distributed on a circle. 


In the year 1935 Johannes G. van der Corput [ ITT] introduced a general construc¬ 
tion principle of a sequence in [0, 1) which resembles exactly the sequence ([Tj) . This 
construction is based on the reflection of the binary digits of non-negative integers n in 
the comma. If n has binary expansion n = no + n±2 + + • • • (which is of course 

finite) with rij € { 0 , 1 } for j = 0,1,2,... , then the nth element y n of the sequence is 
given by 

_ no ni ri2 
Un 2 ' 2 2 ' 2? ' 

The resulting infinite sequence ( y n ) n >o is nowadays called the van der Corput sequence 
in base 2 or the binary (or dyadic) van der Corput sequence. 

Since its introduction in 1935 a multitude of mathematicians studied the van der 
Corput sequence with respect to its distributional properties, but also many variations 
of the van der Corput sequence have been introduced and studied. A search on American 
Mathematical Society’s MathSciNelQ using the keyword “van der Corput sequence” led 
to 93 hits, and presumably there are many more papers dealing with this subject but 
do not explicitly mention van der Corput in the title. The van der Corput sequence in 
base 2 can be seen as the prototype of the most important sequences, even in the multi¬ 
dimensional case, which are nowadays used in modern quasi-Monte Carlo simulation 
algorithms. 

On the occasion of the 125th anniversary of the birth of Johannes G. van der Corput 
(Sept. 4, 1890 - Sept. 25, 1975) we give a survey of many interesting results on the van 
der Corput sequence and we present prominent generalizations thereof. 

The paper is structured as follows: in Section [2] we present results on uniform dis¬ 
tribution modulo 1 of the van der Corput sequence. Generalizations and variants of 
the van der Corput sequence are presented in Section [3j In Section [4] we present some 
multi-dimensional generalizations of the van der Corput sequence. 

Before we proceed, we introduce some notation. By N we denote the set of positive 
integers {1, 2,3,...} and we write N 0 = FfU{0}. The natural logarithm is denoted by log. 
By {a:} we denote the fractional part of a real x and ||x|| denotes the distance of x to the 

1 www.eims.org/mathscinet/, last checked on 12th June 2015 
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nearest integer, i.e., ||x|| = min({x}, 1 —{x}). For 6 e N, 6 > 2, let Zb = {0,1,... ,6—1}. 
Furthermore, let G b be the set of all permutations of Zb. 

For functions f,g:N—t M + , we sometimes write g(N) <C /(IV) or g(N) = 0(f(N )), 
if there exists a constant C > 0 such that g(N) < Cf(N) for all N e N, N > 2. If we 
would like to stress that C depends on some parameter, say p, this will be indicated by 
writing -C p or O p . 


2 Uniform distribution modulo 1 and the van der 
Corput sequence 


According to Weyl |143j a real sequence ( x n ) n >o is uniformly distributed modulo 1 if for 
all 0 < a < b < 1 we have 


lim 

AT —>-00 


#{n € (0,1,..., N - 1} 
N 


{x n } e [a, b)} 


b — a, 


( 2 ) 


i.e., the relative number of points in [a, b) converges to the length of the interval. See 
[25, H5| f° r introductions to the theory of uniform distribution modulo 1. 

Two of the first examples of uniformly distributed sequences are the so-called (na)- 
sequences for irrational a, where x n = na (see [143j ). and the van der Corput sequence. 
Although van der Corput introduced his sequence only with respect to binary digits, we 
introduce a first generalization of his idea to 6-adic digit expansions. 


Definition 1 Let b G N, b > 2. 


• The b-adic radical inverse function is defined as Y b : No -4 > [0,1), 


Y b (n) 


no ni 
b + b 2 + b 3 + 


for n G No with 6-aclic digit expansion n = no + nfb + n 2 & 2 + • • •, where n % e 
(0,1,..., 6 — 1}. In other words, Y b (n) is the reflection of the 6-adic digit expansion 
of n in the comma. 


• The van der Corput sequence in base 6 is defined as Y b := (y n )n> o with y n = Yb(n). 


For example, for 6 = 2 we have 


n 

binary n 
binary Y 2 (n) 
Y 2 (n) 


01234 5 6 7 8 

o. i. io. ii. loo. ioi. no. in. looo. 

0.0 0.1 0.01 0.11 0.001 0.101 0.011 0.111 0.0001 

0 1 1 3 1 5 3 7 J_ 

U 2448 8 8 8 16 


It can be shown by elementary counting methods that the van der Corput sequence 
in base 6 is uniformly distributed modulo 1. We give a sketch of the proof here. 


Proof sketch. For fixed me N and k € (0,1,..., b m — 1} an element y n of the van der 
Corput sequence in base 6 belongs to Jk,m. = [U if and only if n = A(k ) (mod b m ) 
where A(k) is a uniquely determined integer from ( 0 , 1, ..., b m — 1}. Hence exactly one 
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of b m consecutive elements of the van der Corput sequence belongs to Jk, m - This implies 
that 


#{n G {0,1,..., N - 1} : y n e 
and hence 



with 6 G {0,1}, 


lim 

N —>-00 


#{n g {0,l,...,iV — 1} : y n eJ k ,m} 
N 


1 

b rn 


length(J fcjm ). 


( 3 ) 

( 4 ) 


Arbitrary intervals [a, b) C [0,1) are approximated from the interior and the exterior 
by finite unions of intervals of the form ,/ fc m . In this way the result (J4]) carries over to 
the general case. This means that the van der Corput sequence in base b is uniformly 
distributed modulo 1. □ 


For a detailed version of this proof we refer to m Proof of Proposition 2.10]. 


2.1 Discrepancy 

Quantitative versions of (|2]) are usually stated in terms of discrepancy. For a real 
sequence X = ( x n ) n > 0 the local discrepancy is defined, for x G [0,1], as 

a _ #{ n 6 {0,1,..., N- !} : {^n}e[0,x)} 

jy X. 

Sometimes we will abbreviate the counting part in the local discrepancy by A(x; N', X) := 
#{n G {0,1} : {x n } G [0,a:)}. 

Then the extreme discrepancy is defined as 

D n (X) = sup \A x ,n{x') - A x ,n{x)\. 

0<x<x' <1 

The star (extreme) discrepancy is defined as 


D* n (X) = sup |Ax,jv(^)|- 

0<x<l 

The star discrepancy of a sequence X can be interpreted as the Loo-norm of the local 
discrepancy. In this vein one often also studies the L p -norm of the local discrepancy 
and then speaks of the L p -discrepancy of a sequence X, i.e., 

L PiN (X) = |A. YiJV (a;)| p dx^ for p > 1. 

The diaphony (introduced by Zinterhof) is defined as 

1 / oo j AT—1 

Fn ( X ) = 2 E PA E exp(27Tima; n ) 

iv y m= 1 n= 0 

The diaphony is linked to the L 2 -discrepancy by the formula of Koksma: 

Ll N (X) = ^ Q - xX + (5) 

The notions of uniform distribution modulo 1 and discrepancy can be extended to 
the multi-variate case in the obvious way. See, for example, pi 12311251 [g5i m mm rioej. 
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Remark 1 (Some general facts on discrepancy) A real sequence X is uniformly 
distributed modulo 1 if and only if limjv-K» (X ) = 0. For any real sequence X and 
any N we have D* N ( X) < D N (X) < 2D* N (X). Furthermore, for any 1 < p < q < oo 
we have L pN {X) < L qiN (X) < D* N (X). It was first shown by W.M. Schmidt [ I281J in 
1972 (see also jSlEB]) that there exists a constant c > 0 such that for every infinite real 
sequence X we have 


ND^(X) > clogIV for infinitely many N G N. (6) 


In a recent paper Larcher [88] showed that one can choose c = 0.0646363 ... which is the 
largest value for c known so far (instead of 0.06011... in [6], resulting from the lower 
bound 0.12022... for D). The star discrepancy also appears in the Koksma-Hlawka 
inequality for estimating the absolute error of a quasi-Monte Carlo algorithm: for all 
real functions defined on the unit interval with bounded total variation V(f), and all 
sequences X in [0,1) with star discrepancy D* N (X) we have 



i N -1 

77 E fM 

iV n=0 


< v(i)D‘ N (X). 


In the multi-dimensional case V (/) is the variation of / in the sense of Hardy and Krause 
(see [85, [TUB ]). 

It is also well known that for every pG [ 1, oo) there exists a positive number c p with 
the property that for every real sequence X we have 

NL PjN (X) > c p ^JlogN for infinitely many VgN. (7) 

For p G (1, oo) this lower bound was shown by Proinov |l!9j based on results of 
Roth [ 125] and W.M. Schmidt [I3D] for finite point sets in dimension two. Using the 
method of Proinov in conjunction with a result of Halasz [61] for finite point sets in 
dimension two the lower bound follows also for the Li-discrepancy. It follows from a 
paper of Hinrichs and Larcher [73] that one can choose c 2 = 0.0515599... and from a 
paper of Vagharshakyan [ 110] that one can choose C\ = 0.01138... (see also [80] for a 
further discussion). 

Both lower bounds <§ and ([7]) are optimal with respect to the order of magnitude 
in N. 


2.2 The discrepancy in the binary case 

In this section we collect results on the discrepancy of the classical van der Corput 
sequence in base 2. Many of these results are valid in analogous form for arbitrary bases 
b > 2 and even for generalized van der Corput sequences. These generalized van der 
Corput sequences will be introduced in Section 21 and then we will also collect results 
on discrepancy in the general case. 


The star discrepancy. A first estimate of the star discrepancy of the van der Corput 
sequence in base 2 was already provided by van der Corput in his 1935 paper ]l4lj . He 
showed that 


ND* n {Y 2 ) < 


logN 
log 2 


for every N >2. 
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Comparing this estimate with ([6]) we find that the star discrepancy of the van der Corpnt 
sequence in base 2 is optimal with respect to the order of magnitude in N. Later many 
authors worked to improve the discrepancy bound and to find the smallest possible 
constant in the log TV-term. The first who found this optimal constant were Haber jCT] 
and Tijdeman (unpublished according to [%K]). Here we state the slightly more explicit 
result of Tijdeman. 

Theorem 1 (Tijdeman) For every N G N we have 

ND* n (Y 2 ) < 1< ^ V + 1. (8) 

NK J ~ 3 log 2 v ; 

Further very exact results were shown by Bejian and Faure [Bj. Their main finding 
is the formula 


nd* n (y 2 ) = nd n (y 2 ) = 5 ; 

r =1 


N 

2 r 


E 


N 

2 r 


N 

2 m 


whenever 1 < N < 2 m . 


(9) 


This formula gives rise to a very simple generating recursion for D* N {Y 2 ). The simpli¬ 
city of this recursion becomes especially apparent when we use the notation D(N) : = 
ND* n {Y 2 ). Then we have 0(1) = 1 and for any iVeN 


D(2N) = D(N), and D(2N + 1) = \(D(N) + D(N + 1) + 1). (10) 


The formula (joj) . or, even better, the formula ( JlO]) , makes it possible to compute D* N (Y 2 ) 
for every iVeN (see Figure |3j) . 



2 4 8 


16 


32 


Figure 3: D(N ) = ND* n (Y 2 ) of the van der Corput sequence in base 2 for N = 
1,2,..., 32, compared with N i->- + 1. Observe how the recursion (10) is reflected 

in the graph, for example, 0(1) = 0(2) = 0(4) = 0(8) = 0(16) = 0(32). 


Based on formula (jo]) Bejian and Faure reproved the discrepancy estimate ([8]). Fur¬ 
thermore they showed: 
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Theorem 2 (Bejian and Faure) We have 


limsup ND* n (Y 2 ) 

N—^O O \ 


log A \ 
3 log 2 J 


4 log 3 
9 + 3 log 2 


and 


1 

M 


M 


E ^(n) 

N=1 


log M 
4 log 2 


+ 0 ( 1 ). 


Using 0. Bejian and Fanre also provided the intervals for which the discrepancy is 
attained: 


Corollary 1 The reals a G [0,1) for which A y 2 ,n((*) = D N (Y 2 ) are those for which 
there exists a sequence (£ r ) r >i of signs e r G {+, —} such that 


a = ~Y J 


1 

2 r 



(mod 1) 


where || jl^ and ||-||(_ are the right and the left derivatives of ||-||, respectively. 

Note that the number of the intervals given in Corollary [l] is the greatest power of 2 
dividing N. 

A further result whose proof is based on (J9| is the following central limit theorem 
for the star discrepancy of Y 2 . 

Theorem 3 (Drmota, Larcher, and Pillichshammer) For every real y and for M —y 

oo we have 


1 

M 


if 


N < M : ND* n (Y 2 ) < 


log N_ 1 I log N \ 
4 log 2 4 V 3 log 2 J 


*%) + o(l), 


where 


®(y) 



exp(—t 2 /2) d t 


denotes the Gaussian cumulative distribution function. That is, the star discrepancy of 
the van der Corput sequence in base 2 satisfies a central limit theorem. 


Proof sketch. It suffices to prove the equivalent limit relation 


1 

M 


if 


N < M : ND* n (Y 2 ) < 


log M y /logM ) 
4 log 2 + 4 V 3 log 2 J 


*%) + o(l). 


Assume hrst that M = 2 m . Using formula ([9]) one can show that 


ND* N (Y 2 ) = Y,W“n + 0(l), 

r =1 


(ii) 


where oj is uniformly distributed on [0,1). Now, by using well known limit theorems 
for lacunary sequences (compare with [5H] and [115], alternatively we can use [ 13(f ) it 
follows that 

m 

s m = Ells'll 

r =1 
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satisfies a central limit theorem. Since EfS'm] = m/4 and Var [/>,„] = m/48 we thus have 

ym /48 

and also all moments converge. Since ND* n {Y 2 ) = S m + 0(1) we get the same limit 
relation for ND* n (Y 2 ) if N is uniformly distributed on {0,1,... ,2 m — 1}. This proves 
the theorem for M = 2 m . The general case, where M is not a power of 2, can be 
reconstructed from the case of M = 2 m by considering moments. A detailed version of 
this proof can be found in [271 Proof of Theorem 2], □ 


Remark 2 It is known from [7J Theoreme 1] that D* N (Y 2 ) = Y 2 ). Hence the presen¬ 

ted results for the star discrepancy of Y 2 remain valid if we replace the star discrepancy 
by the extreme discrepancy. 

Central limit theorems for the star discrepancy of van der Corput sequences in ar¬ 
bitrary bases b can be found in the PhD thesis of Wohlfarter [ 145] . 


The Z.,, -discrepancy. For the //^-discrepancy of the van der Corput sequence in base 
2 the situation is a bit different than for the star discrepancy. The most studied case is 
the case p — 2. Similar to (J9]) Faure [37] proved the following explicit formula: 


4(iVL 2 ^(r 2 )) 2 



N 
2 / 



N 2 
2t 


( 12 ) 


Based on this formula he was able to show the following results for the L 2 -discrepancy 
of the van der Corput sequence in base 2: 


Theorem 4 (Faure) For all N <G N we have 


(NL 2tN (Y 2 )) 2 < 


f log AT 
y 6 log 2 


2 


/ll 2 log3 \ log N 1 
V 3 log2 J 36log2 + 3' 


Moreover 


lim sup (NL 2:N (Y 2 )) 2 

N^OO \ 


/ log N \ 2 (11 2 log 3 \ log A r \ 

^6 log2 ) ~ \Y ~ log2 ) 36 log2 J 

7 11 log 3 / log 3 \ 2 

” 81 + 108 log 2 + y 6 log 2 ) 


A proof of these results can be found in |33, where a similar study is also performed 
for the symmetrized version of the van der Corput sequence (see Section 3.1 for details). 


For p > 1, the following estimates are known. 
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Theorem 5 (Pillichshammer) Let p > 1. For all N G N we have 


c nl PjN (y 2 )Y < 


log N 


+ 0 ((log N) 


p-il 


v 6 log 2 / 

where the constant in the O-notation only depends on p. If N is of the form 

2m+l / / i \ m+l^ 


N = 


we have 


In particular, 


I / 7Yi I 

(NL P , N (Y 2 )Y > - (- + - - (-1) 
NL PtN {Y 2 ) _ 1 


9 ■ 2 r 


lim sup 

N—^OO 


log N 6 log 2 


= 0.2404.. 


(13) 


A proof of these results can be found in | TT7 j . Equation (13) for 1 < p < 2 can also be 
found in [TBJ 1121 ] and for p = 2 in [CEO • 

A comparison of Theorem [5] with (J7]) shows that the L p -discrepancy of the van der 
Corput sequence in base 2 is not of optimal order of magnitude in N. This problem can 
be overcome with a symmetrization technique which will be explained in Section |3 .1 


On average the L p -discrepancy of the van der Corput sequence is large but behaves 
very regularly, as can be seen from the following theorem, which is m Theorem 3]. 


Theorem 6 (Drmota, Larcher, and Pillichshammer) For every p > 1, for every 
real y, and for M —> oo we have 


1 

M 


# 


N<M : NL PtN (Y 2 ) < 


\ogN y j log N \ 
8 log 2 8 V 3 log 2 J 


$ (y) + o(!), 


that is, the L p -discrepancy satisfies a central limit theorem. 


2.3 The discrepancy of the generalized van der Corput se¬ 
quence 

Definition of the generalized van der Corput sequence. This type of sequence 
was introduced by Faure [32] with the aim of improving the constant term in the asymp¬ 
totic behavior of low-discrepancy sequences. 


Definition 2 Let b G N, b > 2, and let E = (cr r ) r > 0 be a sequence of permutations of 
{0,1,..., 6-1}=: Z 6 . 

• The b-adic radical inverse function with respect to E is defined as Y f f : N 0 —» [0,1), 

V E/ v _ ff o(n 0 ) , o-i(ni) , a 2 (n 2 ) | 

h (n) - —— + -gr~ + -gr~ + ■■■, 

for n G No with 6 -adic digit expansion n = no + n\b + n 2 b 2 + • • •, where Uj G Z&. 

• The generalized van der Corput sequence in base b associated with E is defined as 

'■= (yn)n >o with y n = Yfin). 
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If £ = (cr) r >o is constant, we simply write Y b s = Y b . The original van der Corput 
sequence Y b in base 6 is obtained with the identical permutation, i.e., Y b = Y b d . 

It is an easy extension of the classical case of Y b that for every £ the generalized van 
der Corput sequence in base b associated with £ is uniformly distributed modulo 1. See 
also [32, Propriete 3.1.1] for a short proof using a result on 6 -additive functions. 


Functions related to a pair ( 6 , cr). In [32] and [TB] the formulas (J9]) and ( 12 ) were 
fully generalized. To recall these formulas we first need to introduce some definitions 
used in the study of generalized van der Corput sequences. For J C [0,1) and X = 
(x n )n> o let A( J] N] X) be the number of indices n £ {0,1,..., N — 1} such that x n £ J. 
Let G b be the set of all permutations of 7h b . For any a £ & b set 

-70- _ ( y(o) £11) Ylzdil 
^b V b ’ b ’ " ' ’ b ) ' 

For h £ 7L b and x E [(k — l)/6, /c/6), where k £ {1,..., 6 }, define 


Fb,h( x ) 


A([0, h/b ); k] Z b ) — hx if 0 < h < a(k — 1), 

(6 — h)x — A([h/b , 1); /c; Z b ) if a(k — 1) < h < 6 . 


Note that ip% 0 = 0 for any cr £ & b . Further we observe that if 6 = 2, then we only have 
two permutations which give either <^21 = || • || if cr is the identity or LP 21 = — 1 | • || if a 
is the transposition ( 01 ). 

Moreover, the functions ip bh are extended to the reals by periodicity. Based on <p bh 
we define further functions which then appear in the formulas for different notions of 
discrepancy. Put 



max ^ , V’ft’ = max (-^), V’b = V£’ + + , 

0 <h<b ’ 0 <h<b 


?» = !>£»> « = and xl = Hb- bplf- 

h =0 h=0 

Note that ^ = max 0 < h<h , <b \^ hh - ip^ h ,\ and xl = ~ Ylh'f- 


Exact formulas. For a sequence X = (x n ) n >o in [0,1) let 

D%(X) = sup A.y , N (t) and D^(X) = sup (-A x ,N(t))- 

0<t<l 0<t<l 

The following theorems link, with exact formulas, the functions defined above to 
the discrepancy of generalized van der Corput sequences. Note that the infinite series 
in these formulas can in fact be represented as the sum of finite terms and geometric 
series, and hence can be computed exactly (see [32, 3.3.6 Corollairc 1] for D N , D~^, 
analogous proofs hold for L 2 : n and ibv). We first give the formulas for the extreme 
discrepancies obtained in m Theoreme 1]: 

Theorem 7 (Faure) For any N € N we have 

00 / N\ 

ndUy?) = Y.i’V (y ). 
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ND- k (Y?) = ftf-"- (£). 

00 / AA\ 

JVD w (n E ) = E <+' ( ^ ) and 
D’„(Y b -) = max(£>U+), £»(+)) 


Next, we state the formulas for L 2: n,Fn, and L^n obtained in |TB] Theoremes 4.1, 4.2, 
4.3], but in a fixed base 6 instead of in variable bases as in [[TO] : 


Theorem 8 (Chaix and Faure) For any N G N, we have 


3 = 1 


N 


(NL %N (Y?)? = T £ F-' ( v ) + p E (if) (7-), 




N 


N 


6 > 


(NF N (Yf)?= 4 ^Y./ t -'(^) and 

1 00 / At \ 


Asymptotic behavior for the extreme and star discrepancies. Concerning the 
asymptotic behavior of the extreme discrepancy Dn, the following generalization of (J8]) 
to arbitrary bases with fixed permutations was proved in m Theoreme 2]: 

Theorem 9 (Faure) For a 6 Sj and for all N G N we have 


where 


Moreover, 


DnK) < 4 



log N + max 




d(Y h a ) := limsup 

TV —>00 


NDnW) 
log N 


^ 6 ,<7 

log b 


(14) 


Remark 3 In particular, see [321 Theoreme 6], for a = id we have (the best constant 
being obtained for b = 3 with ct 3iid = 1/2 log 3) 


^6, id 


b 2 

4(6 + 1) 
6-1 

4 


if 6 is even, 
if 6 is odd. 


(15) 


It is worth noting that = 0, so that ^ = Wb~ an d hence Dj^fYif) = D* N {Yf) for 
the van der Corput sequence Y&. Equation (15) shows that D]\r(Yb) increases and tends 
to infinity with 6. The question of the existence of a permutation cq, such that d(Y^ h ) 


1 id,- 


is bounded (by an absolute constant) can be answered positively, see Section [33 
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Theorem [9] is the starting point for research on best possible upper bounds for the 


extreme discrepancy Dn, see the end of the present Section [273] In the next theorem, 
we give an analogous result for the star discrepancy D* N in a slightly different form than 
the original version in [321 Theoreme 3]: 


Theorem 10 (Faure) Let Tb G &b be the permutation defined by Tb{k) = b — k — 1 
for all k G and let A be the subset of N 0 defined by A = Uh=i Ah with Ah = 

{H(H — 1),, H 2 — 1}. For any permutation a G &b, let a := Tb° a and let 

((N)r>o • (cr, ey . cr, cx, (j , cr, cr, n*, cr, cr, cr, ex,. . .) 

6 e i 6 e sequence of permutations defined by a r = a if r G A and a r = a if r A. Then 

ND* N (Y b A ) _ a+ ff + cv 


d*{Yfi A ) := lim sup 


TV—>-oo 


log IV 


2 log 6 


( 16 ) 


where 


a 


\ = inf — sup V'^’b’ + and a b — inf — sup V] fiT ( 7 “ ) • 


Remark 4 (On the “swapping” permutation Tb) Let S be a subset of No and let 

cr G &b■ Define the sequence Eg = (oy ) r > 0 by oy = cr if r G S and oy = a = Tb o a if 
r ^ S. Then, see [32, Lemme 4.4.1], we have 


NDMYT", 


■ El' 


J = 1 
J6S 


J=1 

J6S 


(f) 

OO 

i + E C l 

3 = 1 

(N\ 
KbiU 


j£S 


(f) 

OO 

+ E « ,+ 1 

7 = 1 

(N\ 

Kip) 


and 




The permutation Tb swaps the functions ip%’ + and and hence, in order to minimize 
D* n = ma x(D~fi,Df l ), one has to find a set S for which the sums with and 
asymptotically split into two equal parts. This is achieved, for example, by the set 


A = {0,2,3,6,7,8,12,13,14,15,...} introduced in Theorem 10 above. Further, since 
+ Vb’ - = > for an y 5 C N 0 , we have 

dn(yF‘) = dUyT)+Dn(yF) = £ (<+ (fi)+<•- (fi)) = D„(Y b n 

This property is the starting point for research on best possible lower bounds for the 
star discrepancy for the family of Yfi sequences, and for a new generalization of this 


family, see Section 3.2 


Remark 5 For reasonably small 6 , the constants and a h a are not difficult to com¬ 

pute, and for the identical permutation (for which = 0 ) it is even possible to find 


them explicitly: from (15), for arbitrary b we obtain 


d*{Y b A ) = 


. 1 

f b 2 

<id _ J 

i 8(6 + 1 ) log b 

2 log b ) 

i 6-1 

1 

[ 8 log 6 
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if b is even, 


if b is odd. 


(17) 











For b = 2, this result was first obtained by Bejian [2] and then rediscovered in |SU] in 
the context of shifted Hammersley and van der Corput point sets. The result in ( 8 jQ] 
states that 

ND^Y? 1 ) < + cVtofhv, (18) 

for all N G N, where c > 0 is a constant and where the constant 1/(6 log 2) is best 
possible. We thus see that there exists an explicit generalized version of Y 2 which 
reduces the constant in ([ 8 ]) by a factor of 2 , and this reduction is best possible. 

ft is also worth mentioning the following result from [SU] for general sequences £ and 
6 = 2 which reflects the properties of £ in order to obtain a low star discrepancy. For 
m 6 N let 

S m (E) := max (#{0 < r < m — 1 : ay = r 2 }, #{0 < r < m — 1 : ay = id}) 


and 


T m (£) := ff{l < r < m — 1 : oy_i = r 2 and ay = id}. 


The quantities S m (E) and T m (E) can be used to give very precise discrepancy estimates 
for a generalized van der Corput sequence based on £. The following theorem was shown 


m 


Theorem 11 (Kritzer, Larcher, and Pillichshammer) Let £ = (o>) r >o be a se¬ 
quence of permutations of7L 2 , and let Yff be the corresponding generalized van der Cor¬ 
put sequence. Then it is true that 


Sr 


IS) + T, 


,(S) 


48 


4 < max ND’ n (Y?) < 


5 m (S) 2T m (E) 56 


1<N<2™ 


+ 


9 


+ y 


Theorem 11 implies that two properties of £ influence the star discrepancy of the 
shifted van der Corput sequence. On the one hand, it is the number of identical per¬ 
mutations id compared to the number of transpositions r 2 in £, and on the other hand, 
it is also how these are distributed. Obviously, it is favorable to have a situation where, 
for as many m G N as possible, approximately half of the permutations among the 
first m elements of £ is the identity, and simultaneously there should not be too many 

changes between id and r 2 . This observation is perfectly reflected in the star discrepancy 
X]id __ 

2 a presented above in (18). 


estimate for YY A 


Asymptotic behavior of the L 2 -discrepancy and the diaphony. As for Theorem 
[9j Theorcm|10|is the starting point for research on best possible upper bounds for the star 
discrepancy (see the end of Section 2.3). These theorems have analogous formulations, 
resulting from Theorem [ 8 j for the L p -discrepancy and for the diaphony Fy We only 


give an excerpt of a series of results in fixed and variable bases that can be found in [pi), 
Section 4.2]: 

Theorem 12 (Chaix and Faure) For a G &b and for all N e N we have 


. , v<n NL 2 iN (Yf) /V 

t i(K ) := linrsup — 1 • — = , - — 

oo log N 6 log 6 


with /3b )(T = inf sup | — 

j =i 
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and 


/«") := limsup 

iV-s-oc log N 


b 2 log b 


Theorem 12 confirms that the generalized van der Corput sequences Yf do not have 


optimal L 2 -discrepancy with respect to the order of magnitude in N, which would be 


Remark 6 In particular, see Theoremes 4.12 and 4.13], for a = id we have: 


*i(n id ) 


b 2 + b - 2 

8 (b + 1 ) log b 
b 2 - 1 


8 b log b 


if b is even, 
if b is odd, 


(cf. Theorem [4] for b = 2) and 


/(n id ) = 


tt 2 ( 6 3 + 6 2 + 4) 

48(6 + 1) log 6 
tt 2 (6 4 + 26 2 - 3) 

486 2 log 6 


if 6 is even, 
if 6 is odd. 


These exact values for f(pf d ) are about four times smaller than the upper bound ob¬ 
tained by Pages with the von Neumann-Kakutani transformation, see [1131 Proposition 
2.5], 


We now summarize the results from [82], where the L 2 -discrepancy of Yp is con¬ 
sidered. From Theorem [4] we know that 


NL<2 1 n(Y2) ~ 


log N 
6 log 2 ' 


(19) 


Again, the question arises how the L 2 -discrepancy of Yp is influenced by the properties 
of E. Opposed to the star discrepancy, here one may ask not only whether the constant 
in the log TV-term can be improved, but, in view of (Jt]) , also whether one could possibly 
obtain a better order of magnitude. To give precise answers, we proceed in a way 
that resembles the results for the star discrepancy. Let E = (cr r ) r >o be a sequence of 
permutations of Z 2 , and define, for m G N, 


Z m ( E) := ^{0 < r < m — 1 : o> = id}. 

Furthermore, let z m ( E) := Z m ( E) — Using the quantity z m the following average-type 

result was shown in [82] . 

Theorem 13 (Kritzer and Pillichshammer) Let E = (oy) r >o be a sequence of per¬ 
mutations of Z 2 for which A := lim„ woo exists. Then the following two state¬ 

ments are true: 


lim sup 

m—> oo 


l 

2 m 


2 m 


E 

N =1 


(nl 2 , n (yP)) 2 

m 


< 


A 4 
16 + 3' 


15 
















• For any e > 0 we have 


lim 

m—> oo 


#{1 < N < 2 m : NL 2 , n {Y. 2 e ) < (log 7V) 1/2+£ } 


= 1 . 


Considering the result in Theorem 13, one might hope that there is a sequence E 


such that the order of magnitude in N of L 2> jv(Y 2 ) is lower than that of L 2 n(Y 2 ). 
Unfortunately, these hopes are crushed by the following result that was also shown in 


Theorem 14 (Kritzer and Pillichshammer) There exists a positive constant c such 
that for all sequences E = ( cr r ) r > 0 of permutations of 7L 2 we have 

NL 2jN (Y^) > clogIV for infinitely many N e FT. 

As we see from this theorem, we cannot improve the order of the L 2 -discrepancy of Y 2 
by permutations, but we can again ask whether one can at least improve the constants 


in (19). For this question, there is a positive answer which implies that permutations 


never worsen the constant in (19) and can even improve the constant considerably. The 
following result is also due to 


Theorem 15 (Kritzer and Pillichshammer) It is true that 


sup lim sup 

S TV—>00 


NL 2>n (Y* 

log AT 


6 log 2 ’ 


inf lim sup 

^ N— KX> 


NL 2>N (Yj 


< 


log N 20 log 2 ’ 

where the supremum and the infimum, respectively, are extended over all sequences E of 
permutations of Z 2 . 

Regarding the latter result, one can show that the value 2Q ^ 9 is attained for the 
sequence E in which the identical permutation id and the transposition r 2 alternate. It 
is even conjectured that this result is best possible with respect to all choices of E. Note 
that, if the conjecture is correct, this would mark a difference to the star discrepancy, 
where an alternating sequence E has disadvantages due to the role of the quantity T m ( E) 
For the L 2 -discrepancy it seems to be the case that, while S m also plays 


in Theorem 11 


a crucial role, the value of T rn has no influence on the results. 


The state of the art for asymptotic constants related to D N , D* N , L 2tN , and 
F n . Generalizations of van der Corput sequences give the smallest discrepancies and 
diaphony currently known within the class of infinite one-dimensional sequences X in 
[0,1). We need the following quality parameters to present these current “records”. Let 


d := inf(dLY)) = inf (limsup — ^ ) , d* := inf(d*(AN) = inf (limsup 

\ N^oc log N J x y y ” x \ 


NDWcy 

log N 


t] := inf(t 1 (A)) = inf lim sup 

x y N— >oo 


X 


NL 2 , n (X) ' 
log N 


t 2 (X) := lim sup 

N— >oo 


(NL 2<N (X)y 

log AT 
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h = inf [t 2 (X) ) and / : = inf (/ (X )) = inf Aim sup ^ \ ; 

x x x \ TV-s-oo log iv y 

where the inhma are taken over all infinite one-dimensional sequences X in [0,1). 

For the extreme discrepancy Djy, it was shown, in chronological order, that d < 0.38, 
in base 12 by Faure [321 Theoreme 4] in 1981, d < 0.367, in base 36 by Faure (3BJ 
Theoreme 1.2] in 1992, and d < 0.354 in base 84 by Ostromoukhov [ 1121 Theorem 4] in 
2009, in each case for a specific permutation of digits in the specified base. Hence with 
the lower bound from Bejian |6J in 1982, one has the current estimate 

0.12 < d < 0.354. 


For the star discrepancy D^, it was shown, in chronological order, that d* < 0.2236, 
in base 12 by Faure [32., Theoreme 4] in 1981, and d* < 0.2223 in base 60 by Ostromouk¬ 
hov m Theorem 5] in 2009, in each case for a specific permutation of digits in the 
specified base. Hence with the lower bound from Larcher [ 88 ] in 2014 one has the current 
estimate 

0.0646 < d* < 0.2223. 


Concerning the L 2 -discrepancy, the best permutations for the set of Y b sequences, de¬ 
noted by (Jb, satisfy Hindoo ti(Y^ b ) = 0, see [16] Theoreme 4.14 and its proof]. However, 
this result is not surprising since t\ = 0 as can be seen from the study of symmetrized 
sequences (see Section 3.1) where estimates for t 2 will be given. 

Finally, for the diaphony, besides earlier results by Proinov and Grozdanov [mi 
with weak constants (around 15), it was shown, in chronological order, that / < 1.59, 
in base 2 by Proinov and Atanassov [121] in 1988, / < 1.316, in base 19 by Chaix and 
Faure ESI Theoreme 4.16] in 1993 and / < 1.137, in base 57 by Pausinger and Schmid 
[ 1141 Theorem 3.2] in 2010. As to the lower bound, the only result we know is that of 
Proinov ms from the year 1986 which states that / > 0.0002. 


Comparison with (na)-sequences. Concerning the extreme discrepancy of (nay- 
sequences, Ramshaw [124] obtained 


d((na)) 


1 

51og((v/5 + l)/2) 


0.4156... 


with the golden ration a = 1 and it is a widely held belief that inf Q d((na)) = 

d((na)). Concerning the star discrepancy of (na)-sequences, Dupain and Sos [2S]J proved 
that 


inf d*((na)) = d*((ny/2)) 


1 

41og( v / 2 + 1 ) 


0.2836.... 


2.4 The elementary descent method 


Context, statement and comments. In this section, we are going to give an 
illustration—on the simple example of the triadic van der Corput sequence Y :i —of the 


method yielding the results given in Sections |2.3[ |2.6[ |3.1[ and |3.3| for the generalized 
van der Corput sequences Yf 1 . A more involved variant of this method also applies to 
NUT digital (0, l)-sequences in prime base and to a larger family, the so-called NUT 
(0, l)-sequences over Z b , that will be considered in Section 3.2 Before this illustration, 
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it is necessary to state the fundamental lemma resulting from the descent method and 
leading to the results. However, we first need some notation and a discretization lemma 
to prepare this statement. 

For convenience, in the following we define the error (the deviation from uniform 
distribution) as 

E(x,N,X) := NA x ,n(x) and then E([x,y),N,X) := E(y,N,X)~ E(x,N,X), 

where A x,n(x) is the local discrepancy defined in Section [2j Further, we call a rational 
number in [0, 1 ) an m-bit number if it belongs to the set Q (b m ) = {0, -fh ,..., ^A}. 

Lemma 1 (discretization lemma) Let a G [0,1] and let m, N be integers with 1 < 
N < b m . Furthermore, let u, v be m-bit numbers such that u < a < v and v = u + b~ m , 
and let x be the unique element in (1^(0),..., Y E {b m — 1)} such that u < x < v. Then, 
with y(a ) := u if a < x and y{a ) := v if a > x, we have 

E(a, N, Y?) = E(y(a), N , Y E ) + (y(a) - a)N. 


This is a classical discretization property resulting from the definition of y(a) which 
implies the equality of counting functions A(a; N\ Y E ) = A(y(a)\ N\ Y E ). From Lemmajl] 
it follows that it suffices to know the values of the error E(a, N,Y b E ) for a G Q (b m ). 
Exactly these values are given by the following lemma which is the key property in the 
proof of Theorem [TJ 


Lemma 2 (fundamental descent lemma) Let m , N, and A be integers with 1 < 

N < b m and 1 < A < b m , and let A = Ai6 m_1 H-b \ m -ib + A m be the b-adic expansion 

of X. Then 


E 





3 = 1 


( 20 ) 


where the functions £ 3 = £j(X,m, N) are defined inductively by £ m = rj m = X m and, for 
1 < j < m, 


Vj 



where g' means the right derivative of g, 


£j = r)j if 0 < rjj < b , and £j = 0 if rjj = b. 


The proof is a descending recursion for the fe-adic resolution of the argument y = 
A /b m from m to 1. At each step, the 6 -adic resolution of y is decreased by 1 . The differ¬ 
ences between the discrepancy functions in a reduction step is controlled by means of the 
functions ip bh while the relation between the intervals depends on the right derivatives 
of these functions. 

With Lemmas [l] and [ 2 J it is easy to obtain Theorem [Tj For instance, consider Df: 
First, from Lemma [lj we get 

lim sup E(y,N,Yf) = D+(Y l ’ : ). 

m^OO yg Q( b m) 

Then, since ijj b ,+ = max 0 <fe<b (p b h -> f° r an y ' m_ bit number y G Q (b m ) we get from Lemma 
[ 2 ] the upper bound 

m / N\ 

£te,vn s )<X>r‘' + (^). 
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Now, using the algorithm for the construction of the £j's from the A/s (at the end of 
Lemma [ 2 ]) in reversed direction, it is easy to construct a A G {0,1,..., b rn — 1} for which 
the above upper bound is achieved by E(Xb~ m , N,Y^). Finally, letting m tend to 00 
gives the formula for D^- 

Theorem [ 8 ] is not so easy and needs a further lemma ([16] Lemrne 5.3]) to explicitly 
compute the functions £f : this means that recovering the functions <p% and (j)% is more 
involved in this case. 

There are four versions of Lemma[2j The first one was worked out for the study of the 
./^-discrepancy of the original van der Corput sequence in base 2 and its symmetrized 
version m Lennne 2.2]. The second one is the generalization of the preceding one 
to arbitrary variable bases with arbitrary sequences of permutations (the so-called Yp 
sequences) in [10, Lennne 5.2] for the study of both the L p -discrepancy and the diaphony. 
The third and fourth ones were necessary to extend the method to digital (0, l)-sequences 
Lemma 6.2] and a generalization of these sequences that also includes sequences 
Lemma 2]. In Theorems [2] and [?] the descent method is only present in filigree 
since, at that time, the relation between the intervals at each step of the descent was 
not found yet, and it was not necessary to get formulas for the extreme discrepancies. 
For non-French speaking readers who would like to go thoroughly into the details of the 
descent method we recommend [30 Lemma 6.2] where Xfp has to be replaced by Yp 1 
and the formula there by Equation (20). 


Illustration of the descent method. The first b m points of Y^p, 7 / 0 , • • • iUb m can 

be represented on a checkerboard B rn of b 2m squares as black squares according to their 
position b m yk on the abscissa (where [ 0 , 1 ) is identified with [ 0 , 6 m )) and their index k 
on the ordinate for 1 < k < b m . This is possible because Equation (20) deals with m-bit 
numbers X/b m , the error E(X/b m , N,Yp) involves the m-bit part of ijk = Yp(k) only. 
Then, in this framework, the error is b m times the number A of black squares minus the 
area of [0, A] x [0,1V] divided by b m , i.e., 

E( A, N, B m ) := b m E(X/b m , N, Yp) = b m A - XN. 


One step of the descent method can be described, and the idea of (^-functions pointed 
out, by means of this representation as explained in the following special case of Y 3 , for 
the first step (from m to m — 1) and with m = 3, i.e., b m = 3 3 = 27, see Figure [4j 

We distinguish three cases corresponding to 9(1 — 1) < N < 9 i, where i G {1,2,3}, 
and for each of them we distinguish the abscissas A of the form 3 fi (blue lines), 3/i + 1 
(red lines), and 3/i + 2 (green lines). Concerning the three cases for N, the idea is to 
express the error by moving the abscissa A to the nearest multiple of 3, 3/i or 3 (fi + 1), 
say 3u, in such a way that the black squares (27//*,, k) with 9(1 — 1) < k < 91 still 
remain in the new rectangle [0, 3u\ x [0,1V] in the three cases. In this way, the process 
is the same for any N within a given case, and it will give a single formula for the error. 
Further, E(3u, N , B 3 ) = E(N — 9(1 — 1), B 3 ) since E(3u, 9 (i — 1), B 3 ) = 0, and hence 
(cf. Figure [5]) 

E{ 3v, N , B 3 ) = E(v, N - 9(i - 1), B 2 ). 

Now, it remains to compute the difference 6 = E( A, N, B 3 ) — E(3u, N, B 3 ) for each 
of the nine possibilities to make the step from B 3 to B 2 \ 

• If A = 3/i (blue lines) then for any N we have <5 = 0 and u — /i. 
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Figure 4: Checkerboard B 3 corresponding to the first 3 3 elements of Y 3 . 


20 




























































































































































0 


A__ 

0 12 3 4 


5 6 


7 8 


9 


N 


1 

| 2 

3 

4 

5 

6 
7 


9 



Figure 5: Checkerboard B 2 corresponding to the first 3 2 elements of Y 3 . 



Figure 6 : Graph of (blue), (red) and <pg d 2 (green); - 03 d = max((pg d i, (/d^). 
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• If A = 3/i + 1 (red lines) then for 0 < N < 9 we have 5 = 2 N and v — fi + 1, and 
for 9 < N < 27 we have 5 = 27 — N and v = fi. 


• If A = 3/r + 2 (green lines) then for 0 < N < 18 we have 5 = N and v — fi + 1, 
and for 18 < N < 27 we have 5 = 2(27 — N) and v — fi. 


Coming back to the remainder anc 
</? 3 d 2 (x) = min (a:, 2(1 — x)) (see Figure 
lead to the first step of Equation (20), where 1 < i < 3 and e = A mod 3: 


setting <^ 3 d 0 = 0 , (x) = min( 2 x, 1 — x), and 

61) it is easy to check that the three cases above 


E 


'A, 

27’ 


,iv,n E ) =£(g,jv-9(i-.i),n ! 


1 , 

+ E3,£ 


N 

27 


The next steps follow in the same way, but it is necessary to control the intervals 
that appear inductively in the descending recursion which explains the sophisticated 
definition of the £j’s in Equation (20) and the length of the proof of this formula (5 
pages in 


Concluding this illustration for base 3, we give the analog of Theorem [2] for the 
triadic van der Corput sequence Y 3 and the corresponding graph (Figure [T]) showing the 
structure of the discrepancy by means of the function - 03 d , see [3U] . 


Theorem 16 (Faure) For N G N we have 


ND- N (Y 3 ) = nd n (y 3 ) = £ </4 d = E V4 d Er + 5 = 


r= 1 


N 


r =1 


N 


N 


3 r 

at P * (Y ) / l °Z N 3 log 2 
ND ^ Y3 > £ 2tei3 + 4 + 2log3’ 


whenever 1 < N < 3 m , 

( 21 ) 


and 


limsup ND* n (Y 3 ) 

TV—xx) \ 


log N A 
2 log 3 J 


3 log 2 

4 + 2 log 3 


It is even possible to obtain generating formulas like for Y 2 (cf. (10)), but, due to 
the form of the function i/j 3 d , they are not as simple as for Y 2 where = || • || (the 
“hat-function”). Again the recursion becomes particularly apparent when we use the 
notation D(N ) := ND* n (Y 3 ). Set D( 0) = 0, then for any N G N 0 


D(3N) = D(N ), D(3N + \) = \D{N) + \D{N + |) + |, 

D(3N + l) = ± D (N) + lD(N + l) + % , D{3N + |) = D(N + A) + |, 

D(3N + 2) = lD(N + 1) + \D{N + |) + |, 

D(3N + |) = \D{N + 1) + \D(N + i) + |. 

For larger bases, the complexity of such generating formulas increases and reduces 
the interest in searching for them. 
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Figure 7: D(N) = ND* n (Ys) of the van der Corput sequence in base 3 for N = 
1,2,..., 81, in comparison with N ^s|| + | + AgA. 


2.5 The von Neumann-Kakutani transformation 


The van der Corput sequence can also be obtained from the von Neumann-Kakutani 
transformation. This was first observed by Lambert [85i, 187]. 

The b-adic von Neumann-Kakutani transformation T b : [0,1) —> [0,1) is a piecewise 
translation map given by 


T b {x) = x 


1 1 

b k b k+1 


for x G [l - A, 1 - with k = 0,1, 2,.... 


A plot of T b for b = 2 is presented in Figure [ 8 } The von Neumann-Kakutani transforma¬ 
tion is an ergodic and measure preserving transformation T b : [0,1) —y [0,1), that is, for 
every Borel set A, T b A = A only if A (A) = 0 or A(A) = 1 , and A (T b A) = A(A), where A 
denotes the Lebesgue measure. 

Likewise, T b can also be constructed inductively with a so-called cutting and stacking 
construction which is often used in ergodic theory: at first cut [ 0 , 1 ) into b intervals 
Ij 1 '* = [|, Add) for j — 0,1,..., b — 1. Define the transformation T bA : [0, —> [|, 1) as 

the translation of /onto / for j = 0,1,..., b— 1. In the second step cut all intervals 
Ij' 1 into b subintervals of the form /j 21 = [^, J j 1 ) for j = 0,1,..., b 2 — 1. Then the 
transformation T b2 : [ 0 , —> [^, 1) is the extension of T bt i which translates 

onto I b 2 _ b+ j +1 for j = 0,1,..., b— 1. In this way one can inductively construct a sequence 
of mappings T b i,T b ^ } T b3 ,.... The first two steps for b = 2 are illustrated in Figure [9j 
Finally the von Neumann-Kakutani transformation is defined as T b := lim ?woc T b n . See 
also [ITT] . 


Define the iterates of T b by 

T b x = T b T b ~ l x for neN, 
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where we formally set Tfx = x. Then it is easily seen that the orbit of the iterates of 
Tb with starting point 0 is exactly the van der Corput sequence. For example for b = 2 
we obtain (cf. Figure [ 8 ]) 

0 —>■ 1/2 —>• 1/4 —>• 3/4 —>• 1/8 

It is known that TJ, is not only ergodic, but it is even uniquely ergodic. From this it 
follows that the orbit of the iterates of Tj \ is for every starting point x G [0,1) uniformly 
distributed modulo 1. Hence the sequence 

can be seen as a natural generalization of the van der Corput sequence in base b, which 
is obtained for the choice x = 0. A proof of this fact and more information in this 
direction can be found in the excellent survey pH] by Grabner, Hcllckalek, and Liardet. 

Pages m showed that the orbit of the von Neumann-Kakutani transformation 
with arbitrary starting point x G [0,1) has star discrepancy of order of magnitude 
0((log N)/N) which is optimal according to Schmidt’s lower bound (jbj). 


2.6 Bounded remainder intervals for the van der Corput se¬ 
quence 

It follows from the lower bound (| 6 ]) that for every sequence A" in [0,1) the sequence 
(. ND* n (X )) N£n is unbounded as N —* oo. On the other hand the question arises whether 
there exist t G [0,1) for which NAx,N{t) is bounded in N. If yes, then an interval [0, t) 
with this property is called a bounded remainder interval. Bounded remainder intervals 
or, more generally, bounded remainder sets have been investigated in a series of papers 
by W.M. Schmidt, see, in particular, [ 127, 125]. 

For the van der Corput sequence in base b one can exactly characterize those t G [0,1) 
which yield bounded remainder intervals. This was first done by Hcllekalek [57] . 

Theorem 17 (Hellekalek) An interval [0, t) is a bounded remainder interval for the 
van der Corput sequence in base b if and only if t is a b-adic rational, i.e., t G Q(fe m ) 
for some m G Nq. 


The sufficiency of the condition in Theorem [17] is almost trivial. Writing the interval 
[0 ,ab~ m ) with a G {0,1,..., b m — 1} as the disjoint union \JlZl[kb~ m , (k + 1 )b~ m ) it 
follows from (J3]) that, for suitable 9 k G { 0 , 1 }, 


NAy b ,N {oh 


a—1 

= £ 

k =0 


N 


+ 9i 


N -— < 1. 
b m 


The necessity of the condition needs some investigation for which we refer to |57j . Proofs 
of the result can also be found in [511 55115TJ . Quantitative versions of Theorem 17 are 
given in [23 [MIES]. 

We begin with the van der Corput sequence in base 2 and present the result of EZ]. 
For t G [0,1) with base 2 representation t = t\/2 + T 2 / 2 2 + • • • and m G N define 


ft(m ) := #{j < m : Tj r j+ 1 } and g t (m) 


1 

2 


m— 1 


E l|2“«||. 

U=0 


We have the following result. 
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Theorem 18 (Drmota, Larcher, and Pillichshammer) Let t € [0,1) with an in¬ 
finite base 2 representation. Then for e > 0 we have 

lim # {N < 2 m : (1 - e)g t (m) < Ay 2|J y(t) < (1 + e)gt{m)} _ 1 

m—)-oo 2 m 

The functions fi and g* are intimately related by the inequality 

\{ft{rn) + 1) < g t (m ) < ^(/ t (m) + 1) for all t £ [0,1) and m £ N, 
o Z 

see [23 Lemma 5]. So we obtain the following corollary. 


Corollary 2 Let t £ [0,1) with an infinite base 2 representation, 
have 


Then for e > 0 we 


lim 

m—> oo 


Jf{N < 2 m : (1 —e)f t (m) < Ay 2i jy(t)} 


= 1 . 


Note that Corollary [2] is a strong quantitative version of Theorem 17 (for 6 = 2). It 
follows from [23 Theorem 9] that for every iV £ N we have | Ay 2 ,n if) | < /t(Llog 2 ivj) +4 
for every t £ [0,1). This result even holds for more general digital (0, l)-sequences over 
7j2, the definition of which can be found in Section |T2 


Next, we present another quantitative version of Theorem [17] for generalized van der 
Corput sequences in base 6. We first need some notation: For a £ [0,1) and a real y > 0, 
let u a (y, k) be the number of occurrences of k £ 7L h among the first \_y\ + 1 digits in the 
6-adic expansion of a,and set U a (y) = u a(y, k). Further we call V a (y) the number 
of pairs (0,6 — 1) or (6 — 1, 0) among the first [y\ + 1 digits in the 6-adic expansion of a. 
Let also log b denote the logarithm to base 6, i.e., logger = (logo;)/log6. The following 
result is [ST, Theoreme 1]: 

Theorem 19 (Faure) For any sequence Yfi we have the following estimates. 

1. For every N £ N we have 

N\A Yf , K (a)\ < t±i(/ 0 (log„Af) + ILT^KtlogjJV) + b - +3. 


2. There exist infinitely many N £ N such that 


JV|A rf , w (a)| > if/ 0 (lo gi ,Af) + )u(log t Af) - 1. 


As an immediate consequence, we obtain the generalization of Theorem 17 to gen¬ 
eralized van der Corput sequences in base 6. A simpler version exists for Y 2 , see pd, 
Corollairc 2], The proof heavily relies on the properties of ^-functions introduced in 
Section 2.3 for the study of D N (Yf s ). 

There exists yet another formula for the remainder of van der Corput sequences Y, 
that draws an interesting parallel between these sequences and (no:)-sequences. It was 
proved (see m Theoreme 2]) using the famous property of elementary intervals (see 
Definition [4] in Section 3.2) for the van der Corput sequence. 
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Theorem 20 (Faure) For given a € [0,1) and N <E N, let their b-adic expansions be 
a = Y < ^=o a jb~-’~ 1 and N = YtLoNib 1 . Put I(N,a ) = m i n (^*> af) and let W(N,a ) 

be the number of pairs (—, +) in the sequence of signs o/(a —A^)j> 0 (disregarding zeros). 
Then we have 


oo j 

NA YbtN (a) = - E E -W' J ' + W a) + W(N, a). 

j= o*=o 

The structure of this formula for Yj, is exactly the same as that for (na)-sequences 
obtained in jBU| Equation 2.5]. For the sake of completeness, we recall this formula: 

For given x 6 [0,1) and N e N, let x = YfLi x j m j and N = Y?Li UiTh (where y % — 0 for 
sufficiently large i) be their canonical expansions in semi-regular continued fractions (the 
( Pi/qO’s approach a and m* = qia — pi). Put s it j = qimj—p^, I(N,x ) = Yih niin^, yf) 
and let W(N,x ) be the number of pairs (—, +) in the sequence of signs of (xi — yi)i >o 
(disregarding zeros). Then we have 


NA(na) ,N (o) 


oo j 

-EE + /(AT, z) + fF(7V, x). 

j=0i=0 


These formulas underline the relation between van der Corput sequences and (na)- 
sequences through systems of numeration that permit the precise study of their irreg¬ 
ularities of distribution. The simplicity of 6-adic systems underlying van der Corput 
sequences as compared to continued fraction systems for (na)-sequences, appears to be 
the reason for a better understanding, and also a better asymptotic behavior of the 
former. In the same vein, these two families of sequences are included in a much lar¬ 
ger family of sequences of minimal discrepancy, the so-called “self similar sequences” 
introduced by Borel in na. 


From Theorem 20 it is possible to obtain an improvement of m Theoreme 1 and 
Corollairc 2] for the van der Corput sequence Y 2 , see [33J Corollaire 3 and Section 4.5] . 


The paper |3S] is a variant of Bohl’s lemma for (net:)-sequences applied to generalized 
van der Corput sequences. Bohl’s lemma for (na)-sequences asserts that the remainder 
on an interval J C [0,1) only depends the length of J. Its very short proof is based on 
the density of (not) and on the fact that the counting function does not change under 
translations of J by a (mod 1). In contrast to systems of numeration, the variant for 
E sequences is much more involved (see |35J Theoreme 1]), but it leads to the same 
property for bounded remainder intervals, i.e., for any Yff sequence, J is a bounded 
remainder interval if and only if lcngth( J) is a 6-adic rational. 


2.7 Functions with bounded remainder 


The question dealt with in Section 2.6 can be generalized. Let / : [0,1] —y M be a 
Riemann-integrable function with Jq f(x) da; = 0 and let X = (x n ) n >o be a sequence in 
[0,1). Under which condition on / do we have 


N -1 


E f( x ™) < 1? 

71=0 
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The classical example studied in Section 2.6 corresponds to f(x) = l[ 0i t)(x) — t for 

t e [0. l). 

In [68] Hcllckalek and Larcher studied this question for sequences induced by the 
von Neumann-Kakutani transformation T, \ (and also sequences induced by irrational 
rotations). For x G [0,1) define 


TV—1 


Gn(x) := V Z(W 


n=0 


Recall that (T fe n (0)) n > 0 is the van der Corput sequence in base b (see Section 2.5). 

Theorem 21 (Hellekalek and Larcher) Let f G C' 1 ([0,1]) with Jq f(x ) dx = 0 and 
Lipschitz continuous first derivative f. Then we have: 

1. If /(0) = /(1), then sup^y |Gjv(x)| < °° f or a U x G [0,1). 

2. If sup N |GV(x)| < °° f or some x G [0,1), then /( 0) = /(l). 

3. If fifi) < /( 0), f/ien — oo < linr inf at-kx) G*at(0) and limsupjv^.^ Gat(0) = oo. 

4■ If fiX) > /( 0), then — oo = lim inf at-kx> Gat(0) and lim supa^oo Gat(0) < oo. 


Recall that the Theorem of Rademacher states that every Lipschitz continuous func¬ 
tion on an interval is continuously differentiable in almost every point of the interval (in 
the sense of Lebesgue measure). Hence for the functions considered in Theorem 21 we 
obtain (by using partial integration) that for h 0 we have 


m = [ /Me— d* = fi- h /; /'Me— d* = fi /"(Ue— d», 

and therefore \ f(h)\ -C \h\~ 2 . 

This motivates the following definition: let T a be the class of all 1-periodic functions 
/ : M —?• R. with an absolutely convergent Fourier series expansion with [/(//)[ -C \h\~ a 
for all h 0 and Jq 1 /(x) dx = 0. In [23] Dick and Pillichshammer proved the following 
result for the classical van der Corput sequence in base b: 


Theorem 22 (Dick and Pillichshammer) Let If = (y n )n >o be the van der Corput 
sequence in base b. For a > 1 we have 


TV-1 

5Z f(y») 


< i 


for all f G J- a . 


n =0 

Furthermore, there exists a function g G such that 


TV-1 


lim sup 

N —^OO 


E g{y n ) 

n =0 


= OO. 


( 22 ) 


More accurate results in this direction can be found in [23J. 
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2.8 Subsequences of the van der Corput sequence 

In recent years also distribution properties of subsequences (or of index transformed 
versions) of the van der Corput sequence have been studied. These properties are 
intimately related to uniform distribution of integer sequences. For m G N a sequence 
(k n ) n >o of integers is uniformly distributed modulo m if 


r #{n G {0,1,... ,N - 1} : k n = j (mod m)} 
hm -—- 

N— »oo N 


1 

m 


for all j G {0,1,..., m — 1}. 


For information on uniform distribution of integer sequences we refer to [851 Chapter 5, 
Section 1], 

Let Y b = (y n ) n >o be the van der Corput sequence in base b and let {k n ) n > 0 be a 
sequence in No- Then it is well known (see, e.g., m Remark 4.4]) that the sequence 
(yk n )n >o is uniformly distributed modulo 1 if and only if the sequence (k n ) n > 0 is uniformly 
distributed modulo b d for all d G N 0 . 

From this insight together with known results on distribution properties of integer 
sequences (see, e.g., [BE] and the references therein) one can obtain the following results 
which are taken from m ■ 

Let Y b = (y n )n >o denote the van der Corput sequence in base b. Then the following 
statements hold true: 


1. The sequence Y b ’ v = ( y U n+v)n>o f° r (w , v ) G NxN 0 is uniformly distributed modulo 
1 if and only if gcd(u, b) — 1. 

2. For every a G FT, a > 2, the sequence (y n a )n> o is not uniformly distributed 
modulo 1. 

3. Let p n , n G N 0 , denote the (n + l)-st prime number, i.e., p 0 — 2,p 1 — 3,p 2 — 5,.... 
Then the sequence (y Pn ) n > o is not uniformly distributed modulo 1. 

4. Let F n , n G No, denote the n-th Fibonacci number, i.e., F 0 — 0, F l — 1, and 
F n = F n _ i + F n —2 for n > 2. Then the sequence b, Flb = (yF n )n> o in base b is 
uniformly distributed modulo 1 if and only if b = 5 k for some k G N. 

5. Let F n be as above. Then, in any base b , the sequence (|/Lio g F„j)n>i is uniformly 
distributed modulo 1. 

6. Let £ G R \ Q or ( = l/dfor some nonzero integer d. Then, in any base b, the 
sequence (?/[<] )n>o is uniformly distributed modulo 1. 

7. For an integer q > 2 and n G No let s q (n ) denote the g-ary sum-of-digits func¬ 
tion. Then, in any base b, the sequence Y^ q = (y Sq (n))n> o is uniformly distributed 
modulo 1. 


In some cases we also have estimates for the star discrepancy of index transformed 
van der Corput sequences which we collect in the following theorem. 

Theorem 23 1. Let ( u,v ) G N x No with gcd(w, b) = 1 and let Y^' v be the van der 

Corput sequence in base b indexed by the arithmetic progression un + v. Then we 
have 

D%(Yn < D' N (Y b ). 
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2. Let b = 5 e and let Y' b Flb be the van der Corput sequence in base b indexed by the 
Fibonacci numbers. Then we have 


D* N (Y*' b ) « b 


Vn' 


3. Let b,q >2 be integers and let Yf q be the van der Corput sequence in base b indexed 
by the q-ary sum-of-digits function s q . Then 

1 «, D- N (Y‘-) loglogJV 


v/BgAT 6 ' 

If q < 14, then the upper bound can be improved to 

DUYp) 1 


Vlog N' 


4■ Let a G (0,1) and let Yff be the van der Corput sequence in base b indexed by 
\n a J. Then 

1 log N 


N° 


N c 


The first point is [75;, Theorem 6.1], the second one is [ 1181. Theorem 2.1]. The third 
and the fourth points are taken from m- 


2.9 Distribution properties of consecutive elements of the van 
der Corput sequence 

One test for the statistical independence of successive pseudorandom numbers is the 
so-called serial test. This test can be viewed as a way of measuring the deviation from 
the property of complete uniform distribution. For a given sequence X = (x n ) n > 0 and 
a fixed dimension s > 2, consider the overlapping tuples 


x n = (. x n , x n+ i,..., x n+s - 1 ) for n G N 0 


and consider the discrepancy of the resulting s-dimensional sequence = ( x n ) n > 0 . 
For a random sequence Z = ( z n ) n > 0 in [0, l) s the law of the iterated logarithm states 
that 


lim sup 

TV—>• OO 


V2ND* n (Z ) 

Vlog log N 


1 A“—a.e., 


where A s is the s-dimensional Lebesgue measure. This serves as a benchmark for the 
order of magnitude of D* N (X( s l). More information on the serial test and further statist¬ 
ical tests for the goodness of fit for the empirical distribution of uniform pseudorandom 
numbers can be found in the book of Niederreiter [1061 Section 7.2]. 


We now return to the van der Corput sequence. Let Yy^ — (y n )n>o where y n = 
(17('n), Yb(n + 1),..., Yb(n + s — 1)). In [QJ Blazekova showed the following result for 


s = 2: 
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Theorem 24 (Blazekova) For N e N we have 

DUY, (2) ) = max Q (l - l) . \ (l - ^) 2 ) + 0(D- n (Y„)). 


In particular, linijv->-oo D* N {Y^) = max(|(l — |), |(1 — |) 2 ) > 0. Thus the van 
der Corput sequence Y b does not pass the serial test for s = 2 and is therefore not a 
pseudorandom sequence. In [SB] Fialova and Strauch showed that every point of 
lies on the line segment described by the formula 


1 

y = x ~ 1 + ik + 


i 

b k +i 


for x G 



i 

6 fc + 1 


with fceNo 


( 2 ) 

and they calculated the limit distribution of the sequence Yf . 

These hirelings have been generalized by Aistlcitner and Hofer [TJ to general s > 2 : 
let T b denote the 5-adic von Neumann-Kakutani transformation. Dehne the map 7 : 
[ 0 , 1 ) ->■ [ 0 , l) s by setting 


7 (t) ■= (t,T b t,T%t,...,TZ~ 1 t) 


and set 

r:=hW : te [0,1)}. 

The Lebesgue measure A on [0,1) induces a measure u on T by setting 

v{A) = \({t : 7 (t) G A}) for ACT. 

Furthermore, v induces a measure on [0, l) s by embedding F into [0, l) s , i.e., 

ia(B) = u(B fl T) for B C [0, l) s . 

Theorem 25 (Aistleitner and Hofer) The limit measure of Y^ is fi. 


Note that the measure fi is a copula on [0, l] s since every distribution function of a 
multi-dimensional sequence is a copula if all uni-variate projections of the sequence are 
uniformly distributed modulo 1 . 

The proof of Theorem 25 in jT] is based on the fact that the van der Corput sequence 
can be viewed as the orbit of the origin under the ergodic von Neumann-Kakutani trans¬ 
formation (see Section 2.5). 


Recently Lertchoosakul and Nair 133 studied the same question for sequences of the 
form 

b(h n T Ui), • • • , Yb(^n T ^s))n> 0 ) 

where (ni,... ,n s ) is a fixed s-tuple of non-negative integers and where the integer 
sequence (k n ) n > 0 is Hartmann uniformly distributed on Z. It is shown in pjfj that the 
asymptotic distribution function exists and is a copula. 


3 Generalizations and variants of the van der Cor¬ 
put sequence 

Several authors have studied different ways of generalizing or modifying the van der 
Corput sequence. In this section, we describe some of the most important variants. 
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3.1 Symmetrized generalized van der Corput sequences 

One way to overcome the defect that, compared to the lower bound in ((7|), the van 
der Corput sequence does not have optimal order of L p -discrepancy is based on a pro¬ 
cess which was initially introduced by Davenport for (na)-sequences (see m or j 2 Hl 
Theorem 1.75]). This method is also known as Davenport's reflection principle. Later 
Proinov [T20] . who named the process symmetrization of a sequence, obtained the same 
result with generalized van der Corput sequences. 


Definition 3 A sequence Y = (y n ) n > o is called symmetric if for all n 6 No, one has 
y 2n + U 2 n+i = 1- A symmetric sequence Y = (y n ) n > o is said to be producedby a sequence 
X = ( x n ) n > 0 if for all n G N 0 one has y 2n = x n or y 2n +i = x n . In this case, the sequence 
Y is called a symmetrized version of X and is denoted as Y = X. 

Two years after Proinov, Faure m Theoreme 3] obtained an exact formula for the 
/^-discrepancy of the symmetrized sequence Y 2 of the van der Corput sequence Y 2 (for 
short, the symmetri zed van der Corput sequence ). The proof is based on the descent 
method (see Section 


2.4) worked out for the first time for that study of Y 2 and Y 2 . 
We remark that the symmetrized van der Corput sequence in base 2 


can also be 

obtained with the help of the tent-transformation T : [0,1] —>- [0,1], given by T(a;) = 
1 — |2x — 11, which is a Lebesgue measure preserving function. If Y 2 = (y n )n> o denotes 
the van der Corput sequence in base 2 , then 

Y 2 = *(D 2 ) = (*(3/„))n>0- 

Theorem 26 (Faure) For any positive integers m and N such that 1 < N < 2 m , we 
have 

m N 2 N 2 

(N L 2)N (Y 2 )) 2 = ^(1 — 2\\2 j Y 2 (N + 1)||) - + 

j = 1 


2 ? 


3 ■ 4 r 


Corollary 3 The preceding theorem implies the following formula and estimate: 

1 


t 2 (Y 2 ) = inf max 

m> 1 l<N<2 m 


m log 2 


E(l-2||2W 2 (Ar+l)||) 


N 

2d 


and 


0.089... = 


421 


6750 log 2 


< t 2 (Y 2 ) < 0.103. 


Numerical computations suggest that the lower bound in Corollary [3] is the exact 
value of t 2 (Y 2 ). 

The proof of the preceding result is based on the descent method which permits a 
very precise study of the /^-discrepancy of Y 2 , but which (until now) is ineffective for 
the Lp-discrepancy. The optimal order of L 2 ^(Y 2 ) was also shown in [9TJ by means of 
a Walsh series expansion of the local discrepancy and Parseval’s identity. On the other 
hand, quite recently, Kritzinger and Pillichshammer [84] obtained a general optimal 
estimation for the L p -discrepancy of Y 2 (with non-explicit constants) by means of strong 
analytic tools. This is a new illustration of (at least) two complementary approaches 


in discrepancy theory that are source a of emulation and progress (see also Sections 2.5 


and 2 . 6 ). 
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Theorem 27 (Kritzinger and Pillichshammer) For every p G [1, oo) we have 


L p , N {Y 2 ) ^p 


V^ogN 

N 


We conjecture that an analogous result to Theorem 27 also holds for 
(this is work in progress). The proof of Theorem 27 is based on dyadic 
and on Littlewood-Paley theory. We sketch the basic ideas: 

A dyadic interval of length 2~i, j G No, in [0,1) is an interval of the 


general bases b 
Haar functions 

form 


/ 


Ij,m 


~m m + 1\ 

—,-:— for m 

Y2F 2i J 


0 , 1 ,.. -, 2 ? 


1 . 


We also define /_ ii0 = [0,1). The left and right halves of I = Ij_ m are the dyadic 
intervals I + = I+ m = Ij+i^m and I~ = I~ m = Ij + i, 2 m +i, respectively. The Haar function 
hi = hj }in with support / is the function on [ 0 , 1 ) which is +1 on the left half of I, —1 
on the right half of / and 0 outside of /. The Loo-normalized Haar system consists of 
all Haar functions hj m with j G No and m = 0, 1,..., 2 J — 1 together with the indicator 
function /x_i 0 of [0,1). Normalized in L 2 ([0,1)) we obtain the orthonormal Haar basis 
of L 2 ([0,1)).' 

The Haar coefficients of a function / G L p ([0,1)) are defined as 


hj,m(f) ■= (/, h j>m ) = j f{t)h jjTn (t) dt for j G N_i and m 

J 0 




where we use the abbreviations N_i := N 0 U { — 1}, % := {0,1,..., 2 J — 1} for j G N 0 , 
and D_! := { 0 }. 

In order to prove Theorem [27] we apply the Littlewood-Paley inequality which involves 
the square function 


S(f) := 


1/2 


E E 

y/GN-i m£Dj 


2 2m ax{0j } 


for / G L p ([0,l)). 


Lemma 3 (Littlewood-Paley inequality) Let p G (l,oo). Then there exist c p ,C p > 
0 such that for every f G L p ([0,1)) we have 

c P \\f\\ Lp <\\S(f)\\ Lp <C p \\f\\ Lp . 

So it remains to find good estimates for the Haar coefficients (/, h j m ) when / is the 
local discrepancy of the first N terms of Y 2 . Let p p fn'' m := (^y 2 N , h j, m .)■ The following 
estimate was shown in [ M] : 


Lemma 4 For j G Nq and m G Oj we have 


\h 


N,sym 

j,m 


if j < 

= 2~ 2j - 2 ifj>\\og 2 N]. 


Furthermore we also have 


\h 


N,sym | 

-1,0 I 



Y 2 (M) < 

N —2 N 


if N = 2 M, 
if N = 2M + 1. 


(23) 
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Equation (23) is the most crucial part in the proof of Theorem 27 It shows how the 
symmetrization trick keeps the Haar coefficient /Efo " 1 small enough in order to achieve 
the optimal L p -discrepancy rate. For the non-symmetrized van der Corput sequence the 
corresponding Haar coefficient is at least of order (log IV)/IV for infinitely many iVeN 
which causes the large L p -discrepancy of Y 2 . 

With the ingredients of Lemmas [3] and [4] one can derive the result in Theorem 27 
(cf. 


Proof of Theorem 27. Using Lemma |3j with / = A~^ N we have 

WU) «„ l|S(A ?2W )|| Ir 


yjSN_i mSBj 


E 2 2 ”"' 0 ’' 1 E (pST) 2 n. 


1/2 


N_i 




1/2 


"p/2 


< 


y 2 2 max foj} 


i ieN-i 


N,sym\2 -< 

(hj,m ) •/ 


mG I 


1/2 


A/ 2 / 


where we used Minkowski’s inequality for the L p / 2 -norm. Now Lemma [4] gives 


^ ^ 22max{0, > 7'} 

iGN_i 


^ V f^j,m ^ ~Ij,r 


raEl 


riog 2 a^i-i 


-p/2 


— 4 _/y 2 II-*-[ 0 > 1 )II- l p / 2 n 2 E 


E i; 


i 

16 


E 2_2j ' 

3= riog 2 N 1 


E !/■ 

Z_/ 1 J,n 


mSD,- 


m£D,' 


"p /2 


"p/2 


14 11 

+ — (l 0 g 2 lV+ 1) + 


'41V 2 IV 2 


12 4riog 2 /VI 




log IV 
IV 2 ’ 


where we regarded the fact that (Emeny, 1 i :j . m = 1 for a fixed j E Nq. 


□ 


A general approach to obtain optimal upper bounds for the L 2 -discrepancy of gener¬ 
alized van der Corput sequences in any base b comes from an inequality between Lx n(X) 
and the diaphony F N (X), see jl22 . Theorem A] and also jTBJ Section 4.2.3]: 

Theorem 28 (Proinov and Grozdanov) There exist two constants cq and a 2 such 
that 

{NL 2)N {X)) 2 < 


( 1 

N 

\7T 

_~2_ 


F[n/ 2 \ (X) + ai J + a 2 . 


Corollary 4 The asymptotic behavior of L 2:N {Y'fp) is related to the asymptotic behavior 
of F N (Yfp) by the inequality 

MU) = lim sup < i f(Yf). 

iV-s-oo log IV 7T 
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The state of the art regarding the optimal constant t, 2 = infx(t2p0) is the following 
(see the end of Section 2.3 for the definition). From Corollary |d| we have the constants 
t 2 = //tt 2 resulting from state of the art for the constants / hence, for short, t 2 < 0.134 
in base 19 by Chaix and Faure in 1993 and t 2 < 0.116 in base 57 by Pausinger and 
Schmid in 2010. However, these values remain larger than the bound obtained by Faure 
in 1990 for base 2 with an exact formula for t 2 (Y 2 ), see Corollary [3J namely t 2 < 0.103. 
This means that 


NL 2N (X ) 

0.051599 ... < inf lim sup — 7’ v 7 < 0.32093 ..., 

X 00 Vlog N 

where the inhmum is extended over all sequences A" in [0,1). Note that the lower and 
upper bounds in the latter equation roughly differ by a factor of 6. 


3.2 (0, l)-sequences and NUT sequences 

General remarks on (0, l)-sequences. We already observed in Section [2] that for 
every m G No and k G {0,1,..., b m — 1} exactly one of b m consecutive elements of the 
van der Corput sequence in base b belongs to the so-called b-adic elementary interval 
[4y, ^44). This fundamental property of the van der Corput sequence is the basis of the 
following definition which goes back to Niederreiter [TDbj; see also [ 24| . 

Definition 4 A sequence (x n ) n > 0 in the unit interval [0,1) is called a (0,1 )-sequence 
in base b if every elementary interval E of the form [r^, (A 1 ), with m G No and k G 
{0,1,..., b m — 1}, contains exactly one element of the point set 

{x n : lb m <n<(l + 1 )b m - 1} for every l G N 0 . 

Of course, for every b > 2 the van der Corput sequence W is a (0, l)-sequence in 
base b. 

To take into account important special cases, it is convenient to work with a slightly 
modified definition. To this end define for every x G [0,1] and for every m G N the 
m-truncated expansion of x by [x]y m = Y%Li£ib~ l , where x = is the 6-adic 

expansion of x, with the possibility that £* = b — 1 for all but finitely many i. Then 
a sequence ( x n ) n > 0 in the unit interval [0,1) is called a (0,1 )-sequence in base b in the 
broad sense if every elementary interval E of the form [A-, OT), with m G No and 
k G {0,1,..., b m — 1}, contains exactly one element of the point set 

{[%n]b,m '■ lb m < n < (l + 1 )b m — 1} for every l G N 0 . 

Faure [321 Proposition 3.1] showed that every generalized van der Corput sequence in 
base b is a (0, l)-sequence in base b in the broad sense. The truncation is required for 
sequences T 6 S where oy(0) = b — 1 for all sufficiently large r. 


ft is known that among all (0, l)-sequences in base b the classical van der Corput 
sequence in base b has the worst star discrepancy. This was first shown by Pillichsham- 
mer [ 117 ] for 6 = 2 and later generalized by Kritzer [77] to general 6 and next by Faure 
[12] to (0, l)-sequences in the broad sense (a further extension to ( t , l)-sequences in the 
broad sense is given in [493, see Section |4.3 for the meaning of t ). 
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Theorem 29 (Faure, Kritzer, and Pillichshammer) Let X b be an arbitrary (0,1)- 
sequence in base b in the broad sense and let Y, b be the van der Corput sequence in base 
b. Then we have 


D’ N (X„) < D' n (Y„). 


In particular, 


ND* N (X b 

hm sup sup — -—- 

iV-s-oo X b log N 


C^b, id 

log U 


where a b>i d is defined as in (15) arid where the supremum is extended over all (0,1)- 
sequences X b in base b. 


For the construction of a (0, l)-sequence in base b one frequently uses the digital 
method which is based on the action of infinite NxN matrices, over a finite commutative 
ring with identity and with b elements, on the 5-adic digits of n e No. For simplicity, 
here we only present the case where b is a prime power and where R = ¥ b , the finite 
held with b elements. For the general case (b > 2 and R an arbitrary commutative ring 
with identity and with b elements) we refer to [ 104L IPS ]. We also require a bijection 
: 7L b —> Ft with the property that 99 ( 0 ) = 0. 


Definition 5 Let b be a prime power. Choose an NxN matrix C = (c rtk ) r , k >0 with 
entries in ¥ b such that for any m € N every left upper m x m sub-matrix of C is 
non-singular. Then the sequence Xfi = (x n ) n > 0 defined by 




00 


E 


T (£n,r) 
b r +l ’ 


OO 

where £ n>r = c Tik ip(n k ), 

k =0 


(24) 


and where the n k 

Ffc. Note that the second summation in (24) is finite and performed in ¥ b , but the first 
one can be infinite with the possibility that <£> -1 (£ Tl)T .) = b 


G 7L b are the 6-adic digits of n, is called a digital (0, 1)-sequence over 

formed in ¥ b , but the first 
1 for all but finitely many r. 


If b is a prime number, then we identify ¥ b with lL b equipped with addition and 
multiplication modulo b. In this case we choose for ip the identity. We then obtain the 
original van der Corput sequence Y b with the NxN identity matrix /, i.e., we have 

n = xi 

It is well known that a digital (0, l)-sequence over ¥ b is a (0, l)-sequence in base b if 
C is such that for all k > 0 we have c T)k = 0 for r sufficiently large (see [2411951 106j). but 
in any case it is a (0, l)-sequence in base b in the broad sense (see [12])• The truncation is 
required for sequences Xf in the case when the matrix C yields digits Lp^lpXn^) = b — 1 
for all sufficiently large r. 


NUT digital (0, l)-sequences. An important special case of digital (0, l)-sequences 
is obtained when the associated matrix C is a non-singular upper triangular (NUT) 
matrix. In this case the first summation in (24) is finite as well. These sequences are 
called NUT digital (0,1)-sequences over F fe . 


In (3N[ Faure proved formulas in the vein 
of Theorems [7] and [8] for the different notions of discrepancies of NUT digital (0,1)- 
sequences over ¥ b (where b is a prime number). 

To recall these formulas we first need to introduce some definitions to deal with 
NUT digital (0, l)-sequences. Let again b be a prime number, identify F& with Z&, and 
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let Lp = id. The symbol l±) is used to denote the translated (or shifted) permutation 
corresponding to a given permutation a G (5& by an element v G F& in the following 
sense: 

(a l±) v)(i) := a(i) + v (mod b ) for all i G F 6 , 
and for any r G No we introduce the quantity 

OO 

Or(N) ■■= E <*, kNk (mod b), 

k=r -\-1 


where the N^s are the 6-adic digits of N. Note that Nk = 0 for all k > m if 0 < N < b m , 
thus 6 r (N ) = 0 for all r > m—1 in this case. Finally, for r > 0, we define the permutation 
5 r by 6 r (i) ■— c r ^ r i (mod b) for i G ¥ b and for future use we set A := (5 r )r>o- We can 
now formulate the following result. 


Theorem 30 (Faure) Let X b be a NUT digital (0,1 )-sequence over Ft with prime b. 
Then, using the notation introduced in Section 2P3, for any N G N we have 


DU*?) 


—y w 

N Z ^X b 

3 =1 


Sj — i \&0j _1 ( N ),+ 





DJ,(X?)= ttE* 

1 3=i 


b 


- fN 
W 



D' N (X b c ) = nmx(D+(X b c ),Dk(X b c )), 


1 OO 

(NL 2 ,N(X b )) 2 — t E' 
“ 3-1 


i.NF N( XC )f ^±^Q, 

^_iw%_i(at) (N\ t 1 ^ .Si-iwei-ipv) Z '/N 


V ) + b 2 ^ ^ 


Tb 


IP 


The formulas in Theorem 30 for D^ and F N show that digital (0, l)-sequences gen¬ 
erated by NUT matrices having the same diagonal entries have the same extreme dis¬ 
crepancy and the same diaphony. Therefore, for Dn and Fn generalized van der Corput 
sequences and NUT digital (0, l)-sequences have the same behavior, i.e., D^{X b ) = 
Dn(Y a ) and Fjy(X b ) = F^(Y A ). This is a remarkable feature of NUT digital (0,1)- 
sequences that yields the asymptotic behavior of D N (X b ) and F N (X b ), directly using 
the results in Section 12.31 

For Dfj, Df r , D* n , and L 2 ,n, the formulas are similar to those for generalized van 
der Corput sequences, but the situation is more complicated because of the quantity 
9j-i(N) which depends on N, via the fe-adic expansion of N — 1, and on the NUT 
generator matrix C, via its entries strictly above the diagonal. Very few results are 
available, for example Theorem 32 below for b = 2, or loose bounds for L 2) n hi a very 


special case in base b by Grozdanov 
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Generalization of NUT digital (0, l)-sequences. A further generalization of NUT 
digital (0, l)-sequences, in relation with the search for best possible lower bounds for 
the star discrepancy, has recently been given by Faure and Pillichshammer in [55]. 
The resulting family in arbitrary integer bases—named NUT (0, 1)-sequences over Z ft — 
include both, generalized van der Corput sequences , and NUT digital (0, l)-sequences 
Xf. So far, this is the most general family of van der Corput type sequences. 

The idea is to put arbitrary permutations in place of the diagonal entries of the NUT 
matrix C defining a NUT digital (0, l)-sequence: 

Definition 6 For any integer b > 2, let £ = (oy) r >o £ and let C = {df) r >o,k>r+i be 
a strict upper triangular matrix with entries in Z&. Then, for n £ No, the n-th element 
of the sequence X^' G = (x n ) n > 0 is defined by 

°o x °° 

Xn = UTT with X n,r = Vr{n r ) + C r U k ( m °d &), (25) 

r=0 k=r-\- 1 


where the n^s are the 6-adic digits of n. We call Xf' c a NUT (0,1 )-sequence over Z 6 . 

lie generalized van der Corput 
then we have Xf' C = X G ', 


30 


If all entries of C are zero, then we have X b ’ — Y b , 
sequences. If a r = 5 r for all r > 0 with 5 r as in Theorem 
the NUT digital (0, l)-sequence where C' = {Cr) r >o,k>r is obtained from C by putting 
invertible entries c r r £ Z& on the diagonal. 

Every NUT (0, l)-sequence over Z& is a (0, l)-sequence in the broad sense (the proof 
works like for generalized van der Corput sequences in 


It is quite remarkable that Theorem 30 for NUT digital (0, l)-sequences in prime 
base b is still valid for this generalization (the basic trick is to keep bijections of Z& “on 
the diagonal”, which is the only place where computing inverses is required, all other 
operations concerning entries above the diagonal are performed in the ring Zb): 

Theorem 31 (Faure and Pillichshammer) Let X^’ c be a NUT (0,1 )-sequence over 
Z b where b > 2 is an integer. Then the results of Theorem 30 are still valid for X^’ c 
with permutations oy replacing the permutations S r . 


The proof of Theorem 31 is rather complex and requires a lot of prerequisites, but it 
follows closely the proofs of [30, Theorems 1-4]. See [53] for a sketch where differences 
between both proofs are pointed out. 


Some applications of Theorem 31 will be given later. We now turn to some special 
digital (0, l)-sequences over F 2 . 


Special digital (0, l)-sequences over F 2 . In base b = 2 one knows a NUT digital 
(0, l)-sequence over F 2 with smaller star discrepancy than the van der Corput sequence 
in base 2. This sequence was given in |117L Theorem 3]. 


Theorem 32 (Pillichshammer) For the star discrepancy of the NUT digital (0,1)- 
sequence over F 2 generated by the matrix 


Ci : = 


(l 1 1 
0 1 1 
0 0 1 
0 0 0 

v: : : 
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we have 


0.2885... 


1 „ ND' n ( Af 1 ) 

— —- < lim sup---—- 

5 log 2 N^oo log N 


5099 

“ 22528 log 2 


0.3265.... 


The next example from pE| shows that the van der Corput sequence in base 2, is 
not the worst distributed with respect to the extreme discrepancy Djy among digital 
(0, l)-sequences in base 2 (but Y 2 is the worst distributed with respect to Djy, D* n , Fn 
and L 2N among NUT digital (0, l)-sequences in base b , see Theorem |35| below). 


Theorem 33 (Faure) Let C 2 be the matrix 


C 2 


f i o o o --A 

1 1 0 0 ••• 

1 0 1 0 ••• 

1 0 0 1 ••• 


Then, for any N e N we have 


2D n (Y 2 ) 


5 

2 N 


< D n ( Af 2 ) 


< 2 D n (Y 2 ) 


and 


D‘ N (Y 2 ) - Y < D‘ N (X?) < D‘ N (Y 2 ) = D N (Y 2 ). 

Moreover, the sequence X 2 is the worst distributed among all (0,1)-sequences in base 
2 (in the broad sense) with respect to the extreme discrepancy D at. 


The sequence X 2 2 was already considered by Larcher and Pillichshammer |H]| to 
show that there exists a symmetrized version of a digital (0, l)-sequence in base 2, 
namely X 2 2 , that does not have optimal order of /^-discrepancy. It should be re¬ 
marked that X 2 2 = Y 2 . We see here that X 2 2 has about the same star discrepancy as 
Y 2 , but its extreme discrepancy D N is about twice D N (Y 2 ) and it is the worst among 
all (0, l)-sequences in base 2. Thus, this sequence appears to be one with bad dis¬ 
tribution properties (recall that the symmetrized sequence of Y 2 has optimal order of 
L 2 -discrepancy). 

The /^-discrepancy of NUT digital (0, l)-sequences over F 2 is studied in [2?j. 


Theorem 34 (Drmota, Larcher, and Pillichshammer) Let X 2 be a NUT digital 
(0,1 )-sequence over F 2 generated by a NUT matrix C and let Y 2 be the van der Corput 
sequence in base 2. Then we have 

(NL 2>n (X°)) 2 < ( NL 2 , n (Y 2 )) 2 < + O(loglV) (26) 


and 


lim sup sup 

N- s-oo x$ 


NL 2)N (X%) 
log AT 


1 

6 log 2 


(27) 
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where the supremum is extended over all NUT digital (0,1 )-sequences Xf over F 2 . Fur¬ 
thermore, for the matrix C\ defined in Theorem \3f\ there exists a real c > 0 such that 


-^2,7v(Af) — c 


logiV 

N 


for infinitely many JVeN. 


For central limit theorems for the discrepancy of NUT digital (0, l)-sequences over 
F 2 we refer to m Theorems 5 and 6]. 


The worst distributed sequence among NUT (0, l)-sequences Xf'. Due to 

the theory of (^-functions for the study of van der Corput sequences and their gen¬ 
eralizations, it is possible to bound the discrepancies and the diaphony of sequences 
obtained. 

In [ 401 Lemmas 1-3], the following bounds were obtained: 

C < Kl < Tt an d Xb < X&S 


where a is an arbitrary permutation in & b (the first estimate actually stems from [321 
Section 5.5.4]). According to Theorems [30] and [3l| (where the L 2 -discrepancy is expressed 
by means of Koksma’s formula (|5|)), we obtain the following theorem (not published until 
now): 


Theorem 35 (Faure and Pillichshammer) For any N e N we have 


D* N (Xf’ c ) < D n {X?' c ) < D* N (Y b ) = D N (Y b ), 

L 2 ,n( Af’ C ) < L 2>N (Y b ) and F N (Xfi’ c ) < F N (Y b ). 

In other words, the van der Corput sequence in base b is the worst distributed among 
NUT (0, 1)-sequences Xf’ , especially among generalized van der Corput sequences Yfi 
and among NUT digital (0,1 )-sequences Xf, with respect to D* N , Djy, L 2 ,n and F jy. 


Remark 7 In [4T) ; Theorem 2 and the subsequent Remark] Theorem [35] was already 
stated for Yfi and Xfi. The first part of Theorem 
dependently of [40| Theorem 2] by using Walsh series analysis. 


for L 2)N (Af) was obtained in- 


For D* n , the result of 


Theorem [35] generalizes the bound from Theorem [29] This is not true for D ; y, according 
to Theorem [1 


The latter theorem, for base 2, has an analog in base b but with a 
non-digital (0, l)-sequence in the broad sense, see [421 . Theorem 5.3]. 


Some applications of Theorems 30 and 31 We present two applications, the first 


one concerning best possible lower bounds for the star discrepancy and the second one 
concerning linear digit scramblings, which are also relevant for applications in quasi- 
Monte Carlo rules, and scramblings of (digital) NUT (0, l)-sequences. 


• Best possible lower bounds for the star discrepancy are dealt with in [53, Section 
5] for the (currently) largest family of van der Corput type sequences which was 
introduced in Definition [6] We only give an excerpt of the results. 
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Theorem 36 (Faure and Pillichshammer) Let Csut be the class of all strict 
upper triangular N x N matrices, let a e &b such that D* N {Yff) = D^{Y^) and let 
£ = (cr, a,.. .). Then (see Theorem 10 for notation) we have 


inf <f (Xf- C ) = 

>'( {n.rjjO/r}'’ 2 log b 

CgCsut 


A proof can be found in [55| Corollary 1], Besides the identity id, it is not diffi¬ 
cult to find permutations satisfying the condition D* N (Yff ) = DjyiYf 7 ) from The- 
Tliis can be done, for example, by using intricate permutations (see for 
Section 2.3] for the definition of intrication of two permutations). 


orem 
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instance 

Moreover, a systematic computer search performed by Pausinger returned 26, 58, 
340, and 1496 such permutations in bases 6 , 7, 8 , and 9, respectively. Further, 


the best sequence currently known wit 
Ostromoukhov (see the end of Section 
Hence in base b = 60 we have 


i respect to the star discrepancy found by 
60 


2.31), namely , satisfies this condition. 




inf d*(X fM 

£G{<To,T 6 oocr 0 } N 

CgCsut 


) = 0 . 2222 , 


The case of identity in Theorem 36 is of special interest because a&^d is explicitly 
known for any integer b > 2, see Equations (15) and ©■ 


Corollary 5 With the notation of Theorem 36, we obtain 

1 


inf inf 

b >2 

CgCsut 


d’(Xf' c ) = d’(Yt A ) = 


4 log 3 


= 0.2275. 


This result can be seen as an analog for van der Corput sequences and NUT (0,1) 
sequences of the best lower bound for (nct)-sequences from Dupain and Sos 


mentioned on page 17 So, compared to the best (na)-sequence, we obtain a smal¬ 
ler value for d* with the generalized van der Corput sequence in base 3 obtained 
by alternating the permutations id and t 3 . 


Finally we show that the constant 1/(6 log 2), which is best possible for the star 
discrepancy of any generalized van der Corput sequence according to [EDI Co¬ 
rollary 4] (cf. Remark [5]), is even best possible for the star discrepancy of any 
generalized NUT digital sequence in base 2. The following result is (55 4 Corollary 
3]: 


s c 

Corollary 6 Let X 2 ' be a generalized NUT digital sequence in base 2 generated 
by the NUT matrix C, with permutations £ G {id, r- 2 } N ■ Then, we have 


inf d*(Xf’ c 

se{id,T 2 U 

CgCnut 


1 

6 log 2 


0.2404 


where Cnut is the class of all NUT NxN matrices over F;,. 
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• The notion of linear digit scramblings and scramblings of (digital) NUT (0,1)- 
sequences stems from Matousek m in his attempt to classify the very general 
scramblings introduced earlier by Owen. A linear digit scrambling is a permutation 
of the set Z& of the form 


ir(k) = hk + g (mod b ), 


where /i / 0 and g are given in Z& with b prime. The definition also works for 
any base b, provided that the multiplication by h remains a bijection. If <7 = 0, we 
obtain the so-called multipliers (or multiplicative factors) h of preceding papers on 
NUT digital (0, l)-sequences (see [40] 03]). The additive factor g is a translation 
also called digital shift (see [8T]). 

It is quite remarkable that swapping permutations 7y are linear digit scramblings 
for any base, since n(k) = 6 — 1 — k = (b — 1 )k + b — 1 (mod b ). 


In the following we give some examples showing the effect of linear digit scram¬ 
blings with regard to the asymptotic behavior of the star discrepancy: 


- In base 2 T 2 and the digital shift k t —> k + 1 over F 2 are the same permutations. 


In base 3, all 6 permutations are linear digit scramblings, namely id, id + 1, 
id + 2, 2id, 2id + 1 and 2id + 2 = 73 . All of them have the same asymptotic 
behavior with respect to D* N and there is no improvement with regard to 

d*(Y^ A ), see Theorem 10 and Equation ( JlT| ). 

In base 4 there is no improvement with linear digit scramblings, but the 
permutation a = 3(12) + 1 (where (12) denotes the transposition which 
exchanges 1 and 2) produces the sequence Yf with d*(Y£) = 1/(5 log 2), 
which is the same as for the sequence associated with the two-dimensional 
Halton-Zaremba point set from 


^id 

Moreover, Y 4 A in (17) has also the same 
asymptotic constant 1/(5 log 2), which means that the single permutation a 
is also equivalent to the sequence of permutations in base 4. 


In base 233, which is our best prime base for D* N up to base 1301 (see [40]). the 
phenomenon becomes really apparent: with the two linear digit scramblings 
p{k) = 89 k and n(k) = 89 k + 44 we have 


lim sup 

N —^00 


nd'AySs) 

log AT 


469 + 365 
466 log 233 


0.328..., 


ND* n (Y* 3 ) 368 

Inn sup-—— < - 

n-> ocf log N 233 log 233 


0.289 


) 


ND- n (Y%2) _ 368 

hm sup —---— < ———— 

jv-xjo log N 233 log 233 


0.289..., while 


lim sup 

N —^OO 


log AT 


232 

8 log 233 


5.32.... 


Here, we observe even more: adding the digital shift 44 still improves the 
behavior of the sequence even without using the sophisticated sequence E_ 4 . 
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This fact is in accordance with a remark of Matousek in [1011 p. 540]: “In¬ 
troducing additive terms makes the situation much simpler and more reg¬ 
ular.” Such computations with very large bases have been performed for 
some scrambling techniques applied to Halton sequences and Faure sequences 
(which are examples of (0, s) -sequences), see Section |4j 

3.3 Van der Corput sequences with respect to different digital 
expansions 

Van der Corput sequences have not only been introduced with respect to the classical 

6 -adic digit expansion. Here we present two further prominent examples. 


Cantor expansions. Van der Corput sequences with respect to Cantor expansions 
have been studied by Faure [32] and Chaix and Faure PI- We call B = (bj)j >o with 
bo — 1 and integers b 7 > 2 for all i > 2 a Cantor base and we set Bk b 0 - ■ -bk for 
k G N 0 . 

The special case of ordinary 6-adic expansions, b > 2 an integer, is contained if we 
choose b\ — b 2 — ■■■ — b and hence Bk = b k ■ The main difference between H-adic and 
ordinary 6-adic expansions is that in the general case the i- th digit can take values in 
(0,..., b ^|_i — 1}, which may vary for each i and even become arbitrarily large. Each 
integer n possesses a unique finite representation 

n = 7 i 0 b 0 + nib 0 bi + n 2 b 0 b 1 b 2 H-= 

«>o 

with n t G {0,..., b l+ \ — 1} for i G No- We will call this the B-adic expansion or the 
Cantor expansion of n. Additionally, each real number x G [0,1) has a representation 
of the form 

boh b 0 hb 2 b 0 hb 2 b 3 bp B i+1 ’ 

with Xi G (0,..., bi + \ — 1} for i G Nq. 


Definition 7 Let B be a Cantor base as described above and let E = (cr r ) r >o be a 
sequence of corresponding permutations oy of ( 0 , 1 ,..., b r+ 1 — 1 } for r > 0 . 

• The B-adic radical inverse function with respect to E is defined as : N 0 —» [0,1), 


Yg (n) = + 

for n G No with Cantor expansion n 

{0,1,... ,6 i+ 1 - 1}. 


Pi(wi) & 2 (n 2 ) 

b 2 b 3 

= n 0 B 0 + ni-Bi + n 2 B 2 + • • •, where n 7 


G 


• The generalized van der Corput sequence in Cantor base B associated with E is 
defined as Yf := (y n )n> o with y n = Yg(n). 

Theorems [7] and [ 8 ] extend in a natural way to Cantor bases, see [321 Section 3.4.2, 
Theoreme 1 bis] and [TB], Theoremes 4.1 to 4.3] for their respective statements. In the 
following, we first give a general asymptotic estimate related to the bases bj of the 
Cantor base B, and then we deal with a special construction associated with two pairs 
( 6 , a ) and (c, p) which is useful for building new low discrepancy sequences from other 
ones already known. 
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Theorem 37 (Chaix and Faure) For every Cantor base B for which Yf'j= 1 bj -C n, 
and every sequence £ of corresponding permutations we have 

ND* N (Yjf) « log A r . 

If the sequence £ consists only of the identities, then the condition Yff=\ bj -C n is even 
necessary. 


The sufficiency of the condition on B is m Theoreme 4.5] whereas the necessity in 
case of identities is pm Theoreme 4.7]. 

We now define the notion of intrication of two permutations in two (different) bases 
introduced in [321 Section 3.4.3]: 

Definition 8 The intrication of two pairs (6, a) and (c, p) is the pair ( be , a ■ p ) defined 
by a ■ p(l) = ca{h) + p[k) with l — bk + h, 0 < h < b and 0 < k < c (i.e., we consider 
Euclidean division of / by 6 , 0 < l < be). 

Lemma 5 1. The function t/v P ’ + satifies the relation 

'^bc P ’ + i x ) = i } b ,+ i cx ) + ' 0 e’ + ( a O f or an V x e 
and an analogous relation is also valid for ff h( p '~ and 'f’lf. 

2. Let Vf be the sequence defined by B = ( bj)j>o and £ = ( &j)j>o with b J+1 = b and 
Oj = a if j is even, and h , +1 = c and aj = p if j is odd. Then Yff c p = Y]f. 


If the Cantor base B and the sequence of permutations £ are periodic with the same 
period J, we get Yff = Yfff with Bj = n/=o bj an d ^ J = Il/=o a j (with the product • of 
Definition [ 8 ] which is associative). In particular, if B = (1, b, b, b, .. .) and £ = (cr, er,...) 
are constant, we obtain Yf = Yffff for any n G N. 

Definition |8] and Lemma [5] allow to build by a simple algorithm a sequence of per¬ 
mutations cr& and to show the following theorem, hence positively answering a question 
set in Remark [3] (see [38, Theorem 1.1 and Section 3.1]): 


Theorem 38 (Faure) For any base b > 2, there exists a permutation cp, G &b such 


that 


d{Yff b ) = limsup 

AT—>-oo 


ND N (Y b ° b ) 
log N 


< 


1 

log 2 


/3-adic expansions. Van der Corput sequences with respect to /3-expansions were in¬ 
troduced and/or studied by Barat and Grabner [3], Ninomiya [102, 1110] and Steiner [L33j . 
There are several equivalent ways to introduce such van der Corput sequences. Here we 
follow the presentation in [133] . 

Given a real number (3 > 1, the expansion of 1 with respect to /3 is the sequence of 
non-negative integers (%)j> 1 satisfying 


1 — O.Gqci2 • • • 


a l a 2 

J + ]P + '" 


with CLjdj+i ■ ■ ■ -< ai «2 • • • for all j > 2 . 


(28) 


The symbol -< denotes the lexicographic order of words. For x G [0,1) the /3-expansion 
of x is given by 


x = 0.eie 2 • • • 




with ej€ j+1 ... -< aia 2 ... for all 3 > 1 - 
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Definition 9 The elements of the /3-adic van der Corput sequence Yp = (y n )n >o are the 
real numbers in [0,1) with finite /3-expansion, i.e., 

{y n : n G No} = {0.eie 2 ... : €j€j+i ■ ■ ■ -< cqaq • • • for all j > 1 and 

qq + iq + 2 ... = 000 ... for some I > 1 }, 

ordered lexicographically with respect to the (inversed) word ... 6261 , i.e., for y n = 
0.6162 • • • and y m = O.e^e^ ... we have n < m if there exists some k > 1 such that Ck < e' k 
and (j = e'j for all j > k. 

For example, if /3 0 = 1+ 2 V ^ , then the expansion of 1 is 1 = 0.11000.... In this case 
the initial elements of the /3 0 -adic van der Corput sequence are 


Vo = 0 , 2/1 = 0.1, y 2 = 0 . 01,|/3 = 0 . 001,//4 = 0 . 101,//5 = 0.0001, y 6 = 0.1001, 


This special instance of a /3-adic van der Corput sequence can also be introduced 
with the help of the Zeckendorf expansion: let (F n ) n >_i be the sequence of Fibonacci 
numbers, i.e., F_ 1 = F 0 = 1 and F n = F n _\ + F n _ 2 for n G N. Every positive integer n 
can be uniquely written in the form 

m 

n = a j F j where aj G (0,1} and a/a/+i = 0. 
i=o 

This expansion is called the Zeckendorf expansion of n. Then the elements of Yp 0 with 
/3o = 1+ 2 V ^ are given by 

m 2 m 

Vn = Yl ZF+i if n = F n with 0 < r x < r 2 < • • • < r m and r l+x >n + 1 . 

z—1 P0 i —1 


This definition can also be found in m 3 - 

Ninomiya (l()9j showed that under certain conditions on /3 the sequence Yp has 
optimal order of the star discrepancy: A Pisot number is an algebraic integer for which 
all algebraic conjugates have modulus less than 1. Bertrand jH| and K. Schmidt m 
proved that all Pisot numbers are Parry numbers which means that the expansion 
(28) of 1 is finite or eventually periodic, i.e., ctia 2 ... = Gq ... a^OO ... or Gqa 2 ... = 
aq ... ad- p a,d - p +1 • • • Od, respectively. I 11 this case /3 is the dominant root of the so-called 
/3-polynomial x d — cqa ^ -1 — • • • — ad (with ad > 0) or x d — a x x d ~ l — ■■■ — ad — ( x d ~ p — 


a\x' 


d—p—1 


ad-p) (where p is assumed to be minimal), respectively. 


The following theorem is [ 109 . Theorem 3.1]. 


Theorem 39 (Ninomiya) If (3 is a Pisot number with irreducible /3 -polynomial, then 
we have D* N {Yp) -C (logiV)/iV. 


Bounded remainder sets for the /3-adic van der Corput sequence were studied in 
m- Van der Corput sequences with respect to abstract numeration systems were 
introduced and studied by Steiner |134j . Under some assumptions, these sequences also 
have optimal order of the star discrepancy with respect to W.M. Schmidt’s lower bound 

(§. 
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3.4 Polynomial van der Cor put sequences 

A polynomial version of the radical inverse function was introduced by Tezuka H37]. 
Let b be a prime power and let be the finite field of order b. Furthermore let F 6 ((x x )) 
be the field of formal Laurent series 

OO 

g — y a k x~ k with a k G F& and w G 7L with a w ^ 0. 

k=w 


For g G F & ((ic -1 )) let [g\ be the polynomial part of g. 

First we introduce a polynomial version of the radical inverse function. Choose a non¬ 
constant polynomial p G F&[x] and a bijection 0 : 7L b —* ¥ b with 0(0) = 0. Let q G F&[x] 
and expand this polynomial in base p(x) rather than in x. I.e., find tq, ri,... , rh G F&[x] 
such that 

q(x) = r h {x)p(x) h + • • ■ + r\(x)p(x) + r 0 (x). 

Here h = |_deg(g)/deg(p)J and r 3 = [q/p’\ (mod p) for j = 0,..., h. Now the radical 
inverse Y bip : Fj,[x] —* F^aV 1 )) is defined as 


Y b M 


r 0 (x) ^ _ r h (x ) 

p(x) p(x) h+1 


Let p b : ¥ b ((x x )) ->• [0,1) be given as p b {Y^= w a k x k ) = E£L max (i,™) ^ \ a k)b k - 
We identify n G No with 6-adic expansion n — tiq + nib + • • • + n m -ib m ~ l with the 
polynomial n(x) = 0(n o ) + 0(ni)a; + • • • + 0(n m _i)x m_1 G F fe [x] and vice versa. 


Definition 10 The polynomial van der Corput sequence in base p over F& is given by 
^ = (Vn)n>0 with y n = p b (Y bjP (n)). 


It should be mentioned that this polynomial version of the van der Corput sequence 
over Ft is a special univariate instance of a so-called generalized Niederreiter sequence 
over F& (see P33 Proof of Theorem 1]). Generalized Niederreiter sequences are a very 
powerful construction of low discrepancy sequences and frequently used as sample nodes 
in modern quasi-Monte Carlo rules. For their definition we refer to [231 Section 8.1]; see 
also Section [4731 of this article. 

If we choose in Definition 10 the polynomial p as p(x) = x + a with a G Fj, then 
F b p ° ly is a digital (0, l)-sequence over ¥ b , see [ 137] . 

If b is a prime number, we have ¥ b = Z& and we can choose 0 = id. If we take 
now p(x) = x, then we obtain that Y£° ly = Y b , the classical van der Corput sequence i 
base b. 


m 


4 Multi-dimensional generalizations of the van der 
Corput sequence 

The final section of this survey is dedicated to multi-dimensional variants and general¬ 
izations of van der Corput sequences. These comprise generalized or scrambled Halton 
sequences as well as Niederreiter’s (t, s)-sequences which are the bases of modern quasi- 
Monte Carlo rules for numerical integration, even in very high dimensions (dimensions 
in the hundreds are not uncommon, for example in applications from financial mathem¬ 
atics). 
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4.1 Generalized two-dimensional Hammersley point sets 

Strictly speaking, these point sets are not multi-dimensional generalizations of the van 
der Corput sequence. However, they are intimately related to generalized van der Corput 
sequences and there is abundant literature devoted to these point sets. 


Definition 11 (generalized Hammersley point set) Let b > 2 be an integer, let 
Yj" be a generalized van der Corput sequence in base b as defined in Definition [ 2 ] and 
let m G No- Then the generalized two-dimensional Hammersley point set in base 6 , 
consisting of b m points in [0, l ) 2 and associated to E, is defined by 

<„:={(lf(«),^) : {0,1,..., 


In order to match the traditional definition of Hammersley point sets, which consists 
of m-bit numbers (see page 18), the infinite sequence of permutations E is restricted to 
permutations such that oy(0) = 0 for all r > m, for instance E = (do,..., cr m _i, id, id, id,...). 
Hence, the behavior of Hf m only depends on the finite sequence of m permutations 
<r = (do, ..., d m _]). From now on, we only consider such generalized Hammersley point 
sets denoted by H° m - If a = (d,... , a), then we simply write HI m . 

If we choose in the above definition <jj = id for all permutations, then we obtain 
the classical two-dimensional Hammersley point set in base b denoted by Hb, m H^ m . 
Two pictures of classical two-dimensional Hammersley point sets are shown in Figure [TO] 
(notice that Figures [5] and [5] in Section [2~4 are also representations of the two-dimensional 
Hammersley point sets H 33 and H^)- 



Figure 10: Two-dimensional Hammersley point sets: "H 2 ,8 hr base b = 2 with 2 8 = 256 
points (left) and H 3,5 in base 6 = 3 with 3 5 = 243 points (right). 


Now we need a definition of the local discrepancy in dimension 2. For short, for a 
two-dimensional point set V = {a? 0 ; • • •, *jv-i} of N points in [0, l ) 2 we define the local 
discrepancy (i.e., the deviation from ideal distribution) as 


A v (x,y) 


H([ 0 ,q;) x [0,y);V) 
N 


-xy, 


where H([0, x) x [0, y);V) #{n 6 (0,1, ..., IV — 1} : x n G [0, x) x [0, y)}. 

The link between the local discrepancies of Hammersley point sets (in dimension 2) 
and van der Corput sequences (in dimension 1 ) is almost trivial and has already been 
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used in Section 2.4 We state it again in a lemma as it is the fundamental key that 


permits transferring results for van der Corput sequences to results for Hammersley 
point sets [Ml Lemma 3] . 

Lemma 6 Let m,N, and X be positive integers with 1 < N < b m and 1 < A < b m . 
Then, 

a ( X N \ a ( A \ 

A y_cr | .-, .- | = A V S 


b m b m 


y b TN b „ 


Proof. The proof is almost immediate: from the definition of FL bm the counting functions 
H( PL bm ) an d N, Y b E ) attain the same value, and on the other hand we have 

(note that ff'H bm = b m ), so that the local discrepancies are equal. □ 

From Lemmas [2] and [ 6 j it is possible to obtain (see m for the details) the following 
result. 


Theorem 40 (Faure) For any m G N and any cr = (cr 0 , ..., cr m _i) G (5™ we have, 
with some c m G [0,2], 


l D b rn(ni m ) = max 


max V ib h 

l<N<b m f-J 0 
~ ~ J =1 


j- 1 >+ 



max y ib h 

l<N<b™ f-f 0 
_ ~ 3 =1 


i- 1> 


-IN 



Corollary 7 For any m G N and any cr = (a 0 , ■ ■ ■ ,cr m - i) G & b l we have, with some 
Cm 6 [-2,2], 

D‘ b ~(K J < O t *m(« t ,m) + 

For the star discrepancy of the classical Hammersley point set in dimension 2 we 
have the following sharp result according to De Clerck ra (see also [OS] and [22] for the 
case b = 2 ): 


Theorem 41 (De Clerck) For m > 2 we have: 

• if m > 2 is odd, then 

b m Dlm{n b ,m ) = + (4 + l) - 


• if m > 2 is even, then 


b m D* bm (n b:m ) = m 


b 2 


4(6 + 1) 


+ 


/ 5 26 + 3 \ 

l v 4 + 4(6+l) 2 J 


1 

46 m 



26 + 3 \ 

(& + 1 } 2 /' 


Note that m = (log 6 m )/log 6 and b m = ffPL b ,m- Hence De Clerck’s result (together 
with Corollary [T]) shows that the star discrepancy of the two-dimensional (generalized) 
Hammersley point set is of order Dlm(PL b ,m) = 0((log6 m )/6 m ) = 0((\ogN)/N), which is 
optimal with respect to the order of magnitude in N = b m according to W.M. Schmidt’s 
general lower bound on the star discrepancy of two-dimensional point sets (cf. Re¬ 
mark [I]). We mention that De Clerck ra also provides the “bad” intervals for which 
the star discrepancy is attained. Further results for the star discrepancy of generalized 
Hammersley point sets can be found in [ 431 [721, 80] . The permutations cr can help to 
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minimize the constant which is associated to the leading term m in the star discrepancy 
formulas, see [43] where several results are obtained with the help of the swapping per¬ 
mutation r b applied in various situations, getting analogs in arbitrary bases of results 
from mm in base 2. 


For the L p -discrepancies the situation is quite different. We have: 


Theorem 42 (Faure and Pillichshammer) For any integer base b > 2 and any m G 


Nn we have 


b m Li tb m (FLb, m ) - + 2 + 46^ 


and 


(b m L 2 ^(fH bpn )) 2 = m 2 


'b 2 — 1 
12 b 
1 


+ m 
1 


'3b 4 + 106 2 — 13 b 2 - 1 
7206 2 + 126 


1 

2 b m 


8 4 b m 72 b 2m 

Furthermore, for arbitrary p G [1, oo) we have 

{b m L P}b rn(FL b ,rn)) p = m p + O p M p ~ 1 ), 

where the constant in the O-notation only depends on p and b. In particular, 


v h m L pb 

lim - 

m— t-oo 


’■('H b , r 


log b r ‘ 


A proof of this result can be found in 


b 2 - 1 
126 log b 


Theorem 42 implies that the L p - 
discrepancy for finite p is not of optimal order according to the lower bounds of Halasz, 
Roth, and Schmidt for arbitrary IV-element point sets (which is L Pj n(V) \/logN/N 
in dimension 2, cf. Remark [I]). This defect of the classical Hammersley point sets can 
be overcome when we consider generalized Hammersley point sets. There are various 
papers which provide sequences of permutations which lead to the optimal order of 
L 2 -discrepancy (see, for example, [S21IS31 EH ESI ESI E3 ESI EH H44j ). We provide an 
exemplary result from 1223- 

Theorem 43 (Faure and Pillichshammer) Let cr G {id, T b } m , where r b (k ) := 6 — 
k — 1 for k G Zfe, and let l(m ) denote the number of components of u which are equal 
to id. Then we have 


{b m L2, b AK m )) 2 = 


'6 2 -iy 
6 4 — 1 


((m — 2 l(m)) 2 — m) + 


b 2 - 1 
126 


1 - 


26" 


(21 (m) — m) 


+m- 


906 2 


8 46 m 726 2 " 


In particular, the L 2 -discrepancy of H b m becomes minimal if we choose a sequence of 
permutations cr in which exactly 

m/2 - 2 ifb = 2, 

/ min (m) := { m/2 — 1 if 3 < 6 < 6, 
m/2 ifb>7, 
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elements are the identity in the case of even m, and 


Iminijn) • 


(m — 3)/2 if 2 < b < 3, 


(■m — l)/2 if b > 4, 

elements are the identity in the case of odd m. In any case we have 


= " ( 62 ~S 2 + 13) + 0(1) ' 


(29) 


A generalization of these results and a detailed discussion can be found in EH- 


Equation (29) shows that the optimal permutations lead to an L 2 -discrepancy of order 


0(y/m/b m ) = 0(y/\ogN/N) where N = b m and this is optimal according to the gen¬ 
eral lower bound of Roth H25|. It is interesting that the minimal ./^-discrepancy only 
depends on the number of id-permutations and not on their position within cr. Not 
quite exact as the result for L 2 ,n, but with the same essence, is the next result for the 
Lp-discrepancy (with p G [1, oo)) of TL 2m - 

Theorem 44 (Markhasin) Let p G [l,oo). Let cr G {id, r 2 } m and let l(m) be the 
number of components of cr which are equal to id. Then we have 

^ < ^ l °^ N where N = 2 m 

if and only if |2Z(m) — m\ <C P y/m. 

The proof of this result according to Markhasin [921 ilHD] is indirect via optimality of 
the norm of the discrepancy function in Besov spaces with dominating mixed smoothness 
together with embedding theorems between Besov spaces and Triebel-Lizorkin spaces 
which contain L p -spaces as special cases. The main tool there is the computation of 
Haar coefficients of the discrepancy function, which, for the generalized Hammersley 


point sets, were already computed in EH A direct proof of Theorem [44] via Littlewood- 
Paley theory, which is accessible without knowledge of function space theory, can be 
found in |72| . Although Theorem 44 was proved only in base b = 2 a similar result 


should hold for general bases b. This is work in progress. 

We close this section by mentioning that besides permutations also a symmetrization 
of the classical Hammersley point set leads to the optimal order of L p -discrepancy with 
respect to the order of magnitude in N (see, e.g., iza Theorems 2 and 3]). Furthermore, 
also geometric shifts of the Hammersley point set and their influence on L 2 -discrepancy 
were studied by Bilyk HU- 

4.2 The Halton sequence 

We now come to a very natural extension of van der Corput sequences to arbitrary 
dimensions sGN. 


The classical Halton sequence. Let s G N. A natural generalization of the one¬ 
dimensional van der Corput sequence to dimension s is to choose s different bases 
b\, Z) 2 ,..., b s > 2 and to concatenate component-wise the van der Corput sequences with 
bases bi,b 2 ,..., b s , respectively, to a sequence in the s-dimensional unit cube. Such a 
sequence is called a Halton sequence. 
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Definition 12 Let s > 1 and let b = (61, 62, - - -, 6 S ) be a vector of integers greater than 
1. The s-dimensional Halton sequence in bases 61,62, • ■ ■ , 6 S is given as H b = (y n ) n >0 
where 

Vn = ( y bi W, Y b2 (n),..., Y bs (■ n )) for n e N 0 . 

A minimal requirement on the Halton sequence is that it should be uniformly dis¬ 
tributed modulo 1. This depends on the choice of the bases 61, 62, ..., b s , as can be seen 
in the following theorem. 


Theorem 45 The Halton sequence H b is uniformly distributed modulo 1 if and only if 
bi, 6 2 ,..., b s are pairwise coprime. 


The result in Theorem 45 follows by similar arguments as used in the proof of 
uniform distribution modulo 1 of the van der Corput sequence together with the Chinese 
Remainder Theorem. 

From now on we assume that b is a vector of pairwise coprime integers. In this case 
the Halton sequence even has a low discrepancy of order of magnitude 0((log N) s /N). 
This was first shown by Halton [65]. The following result is due to Niederreiter [IDE] 
and can be shown by a tricky induction argument. Further discrepancy estimates for 
the Halton sequence can be found in [231 SI] 76, SB] US • 


Theorem 46 (Niederreiter) If H b is the Halton sequence in pairwise coprime bases 
61, 6 2 ,..., b s , then 


D’ N (H b ) < 


s 

N 


+ *S 


bj ~ 1 
2 log bj 


log N + 


H 1 ) 


for all NeN. 


It was a long standing open question if the discrepancy estimate for Halton sequences 
is best possible with respect to the order of magnitude in N. Until recently one only 
had the general lower bound due to Kuipers and Niederreiter [55] (which is based on 
the result of Roth [125] for finite point sets) which states that there exists some c s > 0 
such that for every sequence A" in [0, l) s we have 

(log N ) s / 2 

D^(X) > c s - -—- for infinitely many (VeN. (30) 

We add that there is a slight but remarkable improvement of the exponent of the log TV- 
term. It follows from a result of Bilyk, Lacey, and Vagharshakyan ra for finite point 
sets that there exists some c s > 0 and r/. s G (0, |) such that for every sequence X in 
[0, l) s we have 


D* n (X) > c. 


(log N) 


s/2+r/s 


N 


for infinitely many iVeN. 


(31) 


For s > 2 the sharp lower bound for the star discrepancy of infinite sequences is still 
unknown. It is conjectured that it is significantly larger than the one given in (31). In 


general, finding a sharp lower bound seems to be a very difficult problem. We also refer 
to for a discussion of these questions. 

In 2015 Levin [OB] obtained a solution of this question for the important sub-class of 
Halton sequences. He could show that the existing upper bounds on the star discrepancy 
of the Halton sequence are optimal with respect to the order of magnitude in N. 
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Theorem 47 (Levin) Let s > 2 and let B : = 

Then 


sup ND* N (H b ) > 

l<jV<2 mm 0 


m s 

8B 


& 1&2 ■ ■ ■ b 3 and m 0 
for all m > B. 


In particular, there exists some c b)S > 0 such that 


[2B\og 2 b 0 \ + 2. 


D N (dlb) T c b,s 


(log AT)* 

N 


for infinitely many iVeN. 


Atanassov’s results. The bound in Theorem [46] has a slightly defective appearance 
when one considers the constant in the leading log A-term. Consider the quantity 


d* s (X) : = limsup N n D * N ^ 
oo (log N) s 


(32) 


for Halton sequences X = If,. From Theorem [46] it follows that 


< n 


bn ~ 1 


iXi 2 log b j 


= : c,. 


This bound c s is large and grows very fast to infinity for growing dimensions s. For 
example, if bi,...,b s are the first s prime numbers, then the Prime Number Theorem 
implies that bj ~ j log j for large j and hence 


c, = 


n 


n 


J log.? 


n 


J log J 


S! 


=1 2 lo S b j 7=1 2 lo s(i lo g j) 7=1 2 log j 


Also the bounds from |3T, 7B, IM1, T03 J have the same disadvantage. In 2004, Atanassov 
[2] was able to overcome this particular disadvantage (see also the recent survey [45], 
where Atanassov’s method is outlined). 


Theorem 48 (Atanassov) If H b is the Halton sequence in pairwise coprime bases 
h, b 2l ..., b s , then 


D* N {H b ) < 



f lV2j l°g N 

V lQ g & i 



( [bj/ 2J logiV 
V lQ g & i 



1 

N' 


This result implies that 


dl(H b ) < 


1 A IA-/2J 

s! fJi log bj ' 


If again bi,...,b s are the first s prime numbers, then this implies that d* s (H b ) -C j .r 
and hence d*(H b ) tends to zero at an exponential rate as s —> oo. More exact, we have 


lim sup 

s—> OO 


log d*{H b ) 

s 


< — log 2. 


(33) 
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Generalized Halton sequences. A disadvantage of the classical Halton sequence is 
that there seem to exist quite many correlations between the components from the higher 
dimensions, see e.g. [Ml P6] . To overcome this problem one often studies generalized 
Halton sequences which are given by y n = (Y^^n), Y^ 2 (n), ..., Y^ s (n)) for 

n £ No, where Ej = ( &j, r ) r >o are sequences of permutations of {0,1,..., bj} for 1 <j< 


remains valid also for generalized Halton sequences. However, it is important to note 
that the permutations can even lead to a drastic reduction of the star discrepancy of H &. 
For example, in [2J Atanassov constructed specific sequences of permutations E 1; ..., E s 
which lead to exactly such a result: 

Let bi,... ,b s be distinct prime numbers. Then integers k±,... ,k s are called “ad¬ 
missible” to these primes, if bj \ kj for 1 < j < s and if for each set of integers a±,... ,a s 
with bj \ dj for 1 < j < s, there exist integers aq,..., a s such that 

S 

kj 3 P b = dj (mod bj) for all 1 < j < s. 

i=1 

Now let b\,... ,b s be distinct prime numbers and let k±,... ,k s be “admissible” to 
these primes. For each j € {1,..., s} define a sequence of permutations Ej = (r j r )r>o 
of 7L h . by Tj }r (x) := xk r - (mod bj). For these so-called scramblings we have (P, The¬ 
orem 2.3]) 


s; see again 


and the references therein. The discrepancy estimate from Theorem 48 


Theorem 49 (Atanassov) Let bi,... ,b s be distinct prime numbers and let k \,..., k s 
be “admissible ” to these primes. For the sequences of permutations constructed above 
we have 


n* ( 


)< 



E!og bj 

3 = 1 


A ^(l + logfcj) ^ (logN)* f (log N) s A 

(b { - 1) log bj N { N J ' 


From this discrepancy bound it follows easily (see, for example, |48j ) that 

log d* s (H^’-^) 

limsup---= —1, 

s—> 00 S log S 


(34) 


and this is a drastic improvement compared to (33), since now d* s (Hff 1, '"’^ s ) tends to 
zero at a superexponential rate as s —» 00. 


Randomized and deterministic scrambled Halton sequences. Another possib¬ 
ility to overcome the problem of correlations between the components from the higher di¬ 
mensions of Halton sequences is to use so-called “random-start” Halton sequences which 
are based on generalizing the Halton sequence in terms of von Neumann-Kakutani trans¬ 
formations, i.e., the sequence ((T£x 1 , TJfx 2 ,..., TJf a x s )) n > 0 with (x\,X 2 , ■ ■ ■, x s ) £ [0, l) s ; 
see mnnsaEsi- 

One could also consider linear digit scramblings: Various types of linear digit scram¬ 
blings as introduced at the end of Section H have been used by many authors with the 
aim to improve the behavior of Halton sequences for multi-dimensional integration based 
on quasi-Monte Carlo methods (note that also Atanassov’s permutations T 3 . r introduced 
above are specific examples of linear digit scramblings). An extensive comparative study 
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of such methods is given in [45] where a new effective way of scrambling is also proposed. 
This method consists of the search for a list of good multipliers by means of the diaphony 
of one-dimensional projections of Halton sequences, coupled with a criterion based on 
two-dimensional projections that selects the best multipliers in the preceding list (the 
diaphony is considered, rather than the discrepancy, because it is easier and faster to 
compute, see m for more details). Several numerical experiments show a regular and 
better behavior of these multipliers in comparison with other ones (including multipliers 
from the modified Halton sequences of Atanassov). The study also compares linear digit 
scramblings with “random-start” Halton sequences and suggests that these sequences 
have almost the same behavior as original Halton sequences. 


Further remarks. Bounded remainder intervals for Halton sequences are well under¬ 
stood and it is known that they are of the form [0, a\bi mi ) x ■ • • x [0, a s b ~ ms ). We refer 
to Theorem 70] and the references therein. Subsequences of the Halton sequence 
are studied in [7011751 51] . Ergodic properties of /Tadic Halton sequences are studied in 

m ■ 

Also the polynomial version of the van der Corput sequence as presented in Sec¬ 
tion [374] can be generalized to polynomial Halton sequences which are a special sub-class 
of (t , s)-sequences then; see [137]. 


The Hammersley point set. A hnite version of the Halton sequence in dimension 
s — 1 is the Hammersley point set in dimension s (cf. Definition 11 for s = 2). 


Definition 13 (Hammersley point set) Let s,N G N, s > 2, and let b = (6i, 62 ,..., b s - 1 ) 

be a vector of integers with bj > 2. The s-dimensional Hammersley point set in bases 
bi, b 2 , ■ ■ ■ , b a - 1 is given as H b ,N = {y 0 , y x ,..., y N -i} where 



According to a well-known relation between the star discrepancy of infinite sequences 
in dimension s — 1 and hnite point sets in dimension s, whose elements have first co¬ 
ordinate of the form n/N (see, e.g., [241 Lemma 3.45] or [ 106 . Lemma 3.7]) it follows 


from Theorem 46 that if 61 , f) 2 , ■ ■ ■, b s _ 1 are pairwise coprime integers, then 


D* N (Hb,N ) <s,6 


(log N) 
N 


S— 1 


See, e.g., m Chapter 3.4] or [1061 Chapter 3] for more information in this direction. 
For generalized Hammersley point sets we refer to [36]. 


4.3 (t, s)-sequences 

In Section |4.2| we generalized the van der Corput sequence to arbitrary dimension by 
choosing different bases, one per coordinate direction, and concatenating component¬ 
wise van der Corput sequences in these bases. We now introduce a generalization to 
arbitrary dimension where the base b is kept fixed for all dimensions. 
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Definition of (i, s)-sequences. In the one-dimensional case (see Definition [I] in Sec¬ 


tion 3.2), we defined a (0, l)-sequence in base b as a sequence of points such that each 


elementary interval of a certain size contains exactly one element of the corresponding 
segment of the sequence. In higher dimensions, the idea of a (t, s)-sequence in base 
b generalizes this concept. However, the condition that then higher-dimensional ele¬ 
mentary intervals all contain exactly one point of the segments of the sequence is, for 
combinatorial reasons, no longer sustainable in general. Therefore, one needs additional 
parameters to specify how many points lie in the elementary intervals. To be more 
precise, we have the following definition that generalizes Definition |4} The general idea 
of this definition is due to Niederreiter ma. but there also were earlier definitions of a 
similar flavor for special cases by Sobol ’ H32| and Faure 


Definition 14 Let b > 2, t > 0, and s > 1 be integers. A sequence ( x n ) n > 0 in [0, l) s 
is called a (f, s)-sequence in base b (in the broad sense) if for all l,m G No with m > t, 
every elementary interval of volume b t ~ m of the form 


oq cq T 1 


X • ■ ■ X 


a s a s + 1 

fads 7 


_b dl ’ b dl 

with aj,dj G No such that 0 < aj < b dj for 1 < j < s, contains exactly 6* elements of 
the point set 

{[x n }b,m ■ lb m <n<(l + 1 )b m — 1} for all l G N 0 . 

Here the truncation operator [-] 6)m is applied component-wise. 


Some remarks on this definition are in order: 


Remark 8 • The original definition of (i, s)-sequences in base b due to Nieder¬ 

reiter [ 1M . [106 ] is the same but with [x^^m replaced by x n . The definition 
presented here goes back to Niederreiter and Xing (see, e.g., [1081. Chapter 8]). 


• It is well known that every (£, s)-sequence in base b is uniformly distributed modulo 
1 (see, e.g., [Mi EES]). 


• If one would like to have greater flexibility in the definition of a (t, s)-sequence, 
it is also possible to allow the parameter t to depend on the parameter m G N 0 
in the above definition. I 11 this case, the value of t can vary with m, and instead 
of a (t, s)-sequence, one then has a (T. s)-sequence with T = (£(m)) m >o- For the 
sake of simplicity, we concentrate on discussing only (t, s)-sequences with a fixed 
t in the present paper. For further information on (T, s)-sequences, we refer the 
interested reader to p0j and |M, Chapters 4 and 5]. 


Note that in Definition 14 lower values of t mean better distribution properties, 
with t — 0 being the optimal case. For this reason t is called the quality parameter 
of the sequence. For the one-dimensional sequences presented in Section |3.2| a 


value of t = 0 can be reached, however for general sequences in [0, l) s this is not 
always achievable. For instance, a (0, s)-sequence in base b can only exist if s < b 
(if b is a prime power, then the “only if” can be replaced by “if and only if”). The 
question which t for which parameter configuration (s, b) is optimal has gained 
much interest. In [ 107 ] it was shown that for the minimal achievable t we have 


f ^ l ~ l °Sb 


'(b - l)s + b + l' 
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where log fc is the logarithm to the base b. For further information in this direction 
we also refer to m- 

As it was the case for (0, l)-sequences, also for (t, s) -sequences in base b there exists 
a digital version generated by N x H matrices over a finite ring R. As in the one¬ 
dimensional case we restrict ourselves to the case where b is a prime power and R is the 
finite held Fb of order b. Indeed, s-dimensional digital sequences over Fj are obtained by 
choosing, for each j G {1,..., s}, an N x N matrix C l(j> = (c^2)r,fc>o over F&. Using the 
construction method outlined in Definition [ 5 J one obtains a sequence X b {J ' = (x^) n > 0 
for each j G {1,..., s}. The s-dimensional sequence then consists of the points 

x n := (4 1} , ■ ■ ■, x for n G N 0 . (35) 


The sequence (a? n ) n > 0 is often called a digital sequence over F&. 

The matrices C' 1 ' 31 also determine the f-value of the corresponding digital sequence. 
Indeed, let C^\m x m) be the left upper m x m submatrix of C^\ If there exists a 

t G No such that for every choice of nonnegative integers m > t and d\,... ,d s G No with 

d\ + • • • + d s — m — t, 

the hrst d\ rows of C^\m x m), together with 

the hrst <i 2 rows of C^ 2 \m x m), together with 

etc., up to 

the hrst d s rows of C^ s \m x m ), 


(these are m — t vectors in F™) is a linearly independent set over F^, then the digital 
sequence (35) is a (t,s)~ sequence in base b which is then called a digital (t, s)-sequence 
over F(, (see ffence, finding a (t, s)-sequence in base b with a low 

value of t boils down to the problem of finding generating matrices CF 1 ),..., whose 
row vectors satisfy the necessary linear independence conditions. 


The dynamics behind digital (t, s)-sequences was studied by Grabner, Hcllekalek, 
and Liardet. We do not go into details here and just refer to the papers !5U and [62]. 


Generalized Niederreiter sequences. A particular kind of digital (t, s)-sequences 
are so-called Niederreiter sequences, whose points are constructed by means of polynomi¬ 
als over finite helds. These sequences offer an explicit way of constructing the generating 
matrices of powerful digital (t, s)-sequences over F^. They were introduced in [105] . and 
their definition was further extended to generalized Niederreiter sequences by Tezuka in 
m- Following the presentation in [23], a generalized Niederreiter sequence is obtained 
as follows. 

Let s G N, let b be a prime power, and let pi,... ,p s be distinct monic irreducible 
polynomials over ¥ b . Let e 3 := deg(pj) for 1 < j < s. Moreover, choose polynomials yj^ k 
for 1 < j < s, i > 1, and 0 < k < ej with the property that the set {y h i,k : 0 < k < ej} 
is linearly independent modulo p 3 over ¥ b . Then consider the expansions 


yj,i,k(x) 

Pj(xy 


T k, r)x r 1 


r —0 
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over the field F b ((x *)) for 1 < j < s, i > 1, and 0 < k < e r We use these expansions to 
define the generating matrices of a digital sequence; for 1 < j < s let = (c^ )i>i,r>o 
be given by 

cf) = cS j \Q + 1, k, r) G ¥ b (36) 

for 1 < j < s, i > 1, and 0 < k < ej, where i — 1 = Qej + k with suitable integers 
Q = Q(j, i ) and k = k(j , i) with 0 < k < e r 


Definition 15 A digital sequence over ¥ h generated by C^ l \ ..., C^ as given by (36) 
is called generalized Niederreiter sequence. 


Every generalized Niederreiter sequence is a digital (t,s )~sequence over ¥ b with 

t = ~ !)> 

j =i 

see, e.g., [2U ES ODII • 

The concept of generalized Niederreiter sequences comprises earlier constructions 
due to Sobol ’ nm Faure [33], and Niederreiter which are nowadays very popular 
as underlying sample points for quasi-Monte Carlo integration rules (cf. [2T1 Il06j ). 

• A Sobol’ sequence is the generalized Niederreiter sequence where 6 = 2, p\(x) = x, 
and, for j = 2,..., s, Pj(x) is the (j — l)-st primitive polynomial in a sequence 
of all primitive polynomials over F 2 arranged according to nondecreasing degrees. 
Further, there are polynomials 0 ,... ,gj , ej -1 with deg (gj,h) = e j — h + 1 such 
that Hj^k = 9j,k f° r all i G N, k G {0,..., ej — 1}, and j G {1,..., s}, see [1581 
Remark 2]. 

• A Faure sequence corresponds to the case where the base 6 is a prime number 
such that b > s, Pj(x) = x — j + 1 for j = 1, 2,..., s, and all y^k = 1, see [1581 
Remark 3]. In particular, every Faure sequence is a (0, s)-sequence. The matrices 
of a Faure sequence can be given explicitly; there exist s pairwise different elements 
of F b, (3i, ...,I3 S , which are determined by the choice of the pj such that the ith 
matrix is given as 



•• •• ••• / 

where the binomial coefficients have to be taken modulo 6. 


• A Niederreiter sequence corresponds to the case where 6 is a prime power and the 
Vj,i,k ar e of the form yj. t ,k(x) = x k gjj(x ) with polynomials g Jtl G F b [x] such that 
ged (gjj,Pj) = 1, see [ 158 . Remark 4], 
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The currently best constructions for digital sequences with respect to the quality 
parameter t are based on methods from algebraic geometry. These constructions were 
developed in a series of papers by Niederreiter and Xing. For an introduction to this 
subject and an overview we refer to the book by Niederreiter and Xing [ 1U81 Chapter 8]. 
A survey can also be found in the book m Chapter 8]. 


Linear digit scramblings have also been used to improve the quality of (0, s)-sequences 
for quasi-Monte Carlo methods, with the same approach as for Halton sequences in 
Section 4.2, see 


The star discrepancy of (t, s)-sequences. Regarding the discrepancy of (f, s)- 
sequences, it is, not surprisingly, more difficult than in the one-dimensional case to 
come up with good discrepancy bounds. Let us start with discussing the star discrep¬ 
ancy, for which there has been much progress over the past years. Regarding lower 
discrepancy bounds, the best results for dimensions s > 2 are the general bounds (30) 
and (31), excepted in the special case of the two-dimensional (0, 2)-sequence of Sobol’ 
(also a Faure sequence), say S 2 , obtained with the van der Corput sequence Y 2 on the 
x-axis and its image by the Pascal matrix modulo 2 on the y- axis. For this sequence, 
Faure and Chaix [46] were able to obtain the exact order (logIV) 2 for ND* n (S 2 ) with 
a tight lower bound. Together with an upper bound from Pillichshammer in [ 116] . we 
have the following estimate: 


Theorem 50 (Faure and Chaix, Pillichshammer) The sequence S 2 satisfies the in¬ 
equalities 

i .. i 

24(log2) 2 - R" (log/V) 2 - 12(log2) 2 

Until now, this is the only case of a (t, s)-sequence in dimension s > 2 for which the 
exact order of star discrepancy is known. Based on thorough computations, Faure and 
Chaix conjectured that the left inequality should be an equality. 

However, for any dimension, there is a metric result of Larcher and Pillichsham¬ 
mer [23 j for general digital sequences over prime fields. 


Theorem 51 (Larcher and Pillichshammer) Let s € N and let b be a prime num¬ 
ber. Then for almost al^ s-tuples (Ci,. .., C s ) e (F* xN ) s of generating matrices the 
corresponding digital sequences have a star discrepancy satisfying 

n* ^ n \ 0°g N Y log log N f ■ u ■+ i at ^ ™ 

D n > c(o, s)-—- for infinitely many A e N 

with some c(b, s) > 0 not depending on N. 


The situation is different for upper discrepancy bounds, where one can use the struc¬ 
ture of (t, .s)-sequences to show strong results which are better than the metric lower 
bound from Theorem [5TJ Upper bounds for the star discrepancy that hold for arbit¬ 
rary (t, s)-sequences are due to Niederreiter [ HR] . which are based on a clever recursion 
argument. The bounds in [ 104 j were slightly refined in [78], and in a series of papers 

2 Here “almost all” is meant with respect to a certain probability measure y s on the class of all 
s-tuples of N x N-matrices over F&; see IS3I- 
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mmm, further improvements were shown. We shortly summarize some of the most 
up-to-date results. Star discrepancy bounds for (t, s)-sequences can be seen from two 
different perspectives. One may either try to optimize the term involving the highest 
power of logiV, including the constants, as this yields the asymptotically best discrep¬ 
ancy bounds. On the other hand, one may try to optimize the bounds globally such 
that for the non-asymptotic setting the bounds might be better than those that are 
optimized with respect to the leading term. 

Regarding the former approach, Faure and Kritzer showed in [47] the currently best 
general star discrepancy bounds for general (t, s)-sequences, where the constant in the 
leading term is the lowest currently known. The following result was shown in 03 using 
a recursion on s and a discrepancy bound for finite point sets known as (t, m, s)-nets. 

Theorem 52 (Faure and Kritzer) Let s > 2 and let S s be a (t, s) -sequence in base 
b. Then it is true that 


D* N {S s )<c Syb b 


t (log N) s 


N 


+ 0 b f 


(log AT)® -1 ' 
TV 


where 


hW^j(£ib) S ifbiseven, 
M{Mb) S if bis odd. 


If we consider the quantity d* as defined in (32) for (t, s) -sequences in base b, then 
we obtain from Theorem [52] that 


<r s (Ss) < c s , b b\ 


It is known that for any s,b G N, b > 2, one can construct a (t b (s), s)-sequence S* in 
base b with t b (s) x 6 s. These constructions are based on the method of Niederreiter and 
Xing (see |1TT8L Chapter 8]). For such sequences we obtain 


lim sup 

s—>• oo 


logics?) 

slogs 


< -1 


and hence d* s (S*) tends to zero at a superexponential rate as s —* oo. This is the same 
rate as for the special generalized Halton sequence given by Atanassov [2], cf. (34) in 


Section 4.2 Remember that the analog result for Halton sequences given in (33) is much 


weaker. A comparison is also discussed in 


The paper 03 contains also a result where all constants in the star discrepancy 
bounds for (t , s)-sequences are made explicit (cf. [47, Theorem 2]). However, the focus 
of |U] lies on optimizing the leading constant c s . Faure and Lemieux showed in [45). [701 
further upper discrepancy bounds which are not optimized regarding c s but which in the 
non-asymptotic setting (i.e., for smaller values of N) yield better discrepancy bounds 
than m Theorem 2]. The bounds of Faure and Lemieux are as follows. 

Theorem 53 (Faure and Lemieux) Let S s be a (t, s)-sequence in base b. For any 
N e N it is true that 


D* n (S s ) < ^ 



log AT 
log b 



1 

N 


it s-l 

H— y 


b_ 

k\ 



log AT 
log b 
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Theorem [53] was shown in [49] by an adaption of Atanassov’s method for Halton 
sequences to (t, s)-sequences. Finally, in [50( Theorem 1], the following improvement of 
Theorem [53] was shown. 


Theorem 54 (Faure and Lemieux) Let S s be a (t, s)-sequence in base b. For any 
N G N it is true that 


D" N (S.) < ([ 


(log IV 



+ 7 bk 






where 7 fe := 2 — (b (mod 2)). 


Numerical experiments in [HI] for different values of N and the other parameters 
show that the bounds in Theorems 53 and 54 may have an advantage over those in 
Theorem [52] in the non-asymptotic setting. 

Finally, we mention that special discrepancy estimates for generalized Niederreiter 
sequences were recently shown by Tezuka in [ 139] and then improved for even bases by 
Faure and Lemieux m Theorem 2], 


Theorem 55 (Tezuka, Faure, and Lemieux) Let S s be a generalized Niederreiter 
sequence in dimension s. Then it is true that 


D* n (S s ) < 


1 1 hi. n>‘ i - 1 / log A" 

.s! N \ 2 ej [ log 6 

1 ^ b ek + 1 fr (b ei - 1 

+ N ho k\ 2e* 


+ J2 e i 
2—1 

/ log A r 
\ log 5 




where the e v 1 < j < s, are the degrees of the polynomials pj needed in the construction 
of generalized Niederreiter sequences as outlined above. 


The bound in Theorem 55, as well as the earlier result of Tezuka in [ 1~5'9] . 


were 


obtained by relating generalized Niederreiter sequences to the new concept of so-called 
(t, e, s)-sequences as introduced in [139] . Indeed, generalized Niederreiter sequences can 
be viewed as (0, e, s)-sequences, where e = (ei,..., e s ), with the ej exactly corresponding 
to the degrees of the polynomials pj. 


The Lp-discrepancy of (f, s)-sequences. Regarding the L p -discrepancy with p G 
[1, cxi), less is known for general (t, s)-sequences than in the case of the star discrepancy. 
In general, it is clear from the monotonicity of the L p -discrepancies in p that the upper 
bounds on the star discrepancy of (t, s)-sequences are also valid for the L p -discrepancy. 
However, in view of a general lower bound due to Proinov [119] for p G (1, oo), which 
states that there exists some c SjP > 0 such that for every infinite sequence X in [0, l) s 
we have 

(log ivW 2 

L p jv(X) > c s>p -—- for infinitely many N G N, 

the question remains whether there exist (t, s)-sequences for which the L p -discrepancy 
matches this lower bound with respect to the order of magnitude in N. This question 
was answered in the affirmative by Dick and Pillichshammer for p = 2 (see [23 EE]) and 
by Dick, Hinrichs, Markhasin, and Pillichshammer [21] for general p recently. 
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Theorem 56 (Dick, Hinrichs, Markhasin, and Pillichshammer) For p G [1, oo) 

and every dimension s one can explicitly construct a digital sequences S s over F 2 for 
which we have 


J p,N 


(S s )« 


s,p 


(log N) s / 2 
N 


The construction mentioned in the above theorem is based on so-called higher order 
sequences which were introduced by Dick DEI 120], and which are especially useful to 
obtain optimal convergence rates for the integration of smooth functions (see also (24] 
Chapter 15]). One particular instance of such a construction are so-called digit-interlaced 
Niederreiter sequences. 


Conclusion. Nowadays, digital (t, s)-sequences, which have their roots in the simple 
digital construction due to J.G. van der Corput, are among the most powerful construc¬ 
tions of sequences which are used as the sample points for quasi-Monte Carlo integration 
rules 

1 N ~ 1 r 

T 7 Z /(*n) ~ f(x)dx. 

N n =0 

The general link between the distribution properties of the underlying sample nodes and 
the absolute integration error of a quasi-Monte Carlo rule is provided by the famous 
Koksma-Hlawka inequality which was already stated in Remark [l] Introductory texts 
on the quasi-Monte Carlo method and digital sequences are the books |24, H6J H"06] and 
the survey [22] , 
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